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SUMMARY 



A linear program for matrix algebra required ir\ a first course in 

multivariate educational statistics was t^cJuire* throug^self -instruction, 
program is to enable graduate students to acquire ^hrougKseit r 

canonJoal S correl°ticn" 

some knowledge of matrix operations, determinants, linear epene 

vector spaces, and the characteristic equation of the matrix. The linear 
program covers these topics. 

Material from standard texts in matrix algebra was ^^"^^Snee 
a linear program. The material was selected in the light of |\ 

in teaching this subject and the reviews of this mater a n 

variate texts. 

A preliminary trial of the program was carried out with twenty-eight 

graduate students from the University of Chicago and reaction 0 f the 
Hawaii. Analysis of errors made by the students and th . 

students to the material was performed for the purposeof ^vising 
program. The revised material was given to another twenty-nine graduate 
students of the Department of Education to read, and necessary revisio 

was made. 
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INTRODUCTION 



The sequence of educational statistics courses at many 
universities includes an introduction to the use of multivari- 
ate methods in educational research. The courses cover mul- 
tiple and partial correlation, canonical correlation, multi- 
variate analysis of variance, and component and factor analysis* 
Since all textbooks and other literature in the field of multi- 
variate analysis make use of matrix algebra, it is impossible 
to teach the courses without assuming some knowledge of this 
subject on the part of the students. Many schools with a theo- 
retically oriented mathematics department have no course in 
matrix algebra suitable for applied workers. Since it has 
seemed undesirable to include a purely mathematical subject in 
the education curriculum, we have had to depend on our students 
to prepare themselves in matrix algebra by their own resources. 
Although a number of conventional texts are available for this 
purpose (Aiken, 1956 ; Browne, 1958 ; Horst, 1963 ; Ayers, 1962 ), 
the students appear to need more guidance in approaching the sub 
ject than can be offered in a textbook. It was therefore pro- 
posed to prepare a linear program, which will lead the student 
through the material with essentially no errors. 

In recent years several studies have been made of the 
effectiveness of linear programs compared with other methods of 
teaching. The results are still inconclusive, but the linear 
programs definitely take less time than other methods to learn 
the same material. Studies by Hughes and McNamara ( 1961 ), 
Porter ( 1961 ), Hough ( 1962 ), and Smith ( 1961 ) are some of 
the examples which confirm this fact. Bhushan ( 1966 ) also con- 
ducted a study of the effectiveness of two methods of teaching 
elementary matrix algebra and did not find any significant dif- 
ference between the two, but the errorless program took signi- 
ficantly less time than the dialectical program. This result 
encourages me to believe that the student will find the pro- • 
gram an efficient way to prepare themselves in matrix algebra. 




J 



METHODS 



The contents of the program were taken from standard texts in matrix 
algebra (Browne, 1958; Horst, 1963; Hohn, 1958; Kemeny, 1957; Murdock, 1957; 
Schwartz, 1961) and were incorporated in a linear program. Professor Bock 
who teaches a course in multivariate analysis at the University of Chi= 
cago and multivariate texts of Anderson (1958) and Rao (1965) guided the 
selection of the material to be included in the program. The contents 
were divided into eighteen units. The program was revised on the basis of 
the experts* suggestions. 



I 




For the evaluation of the material, it was given to the graduate 
students of the department of Education, thirteen from the University of 
Chicago and fifteen from the University of Hawaii, to read. They were asked 
to read each frame, provide the necessary answer and compare it with the 
one given just after the frame but masked x*ith a piece of paper. If the 
answer did not tally with the given answer, they checked it with a red 
pencil and proceeded further. They were also asked to point out if they 
had any difficulty in understanding any frame. At the end they were given 
an achievement test. The mean achievement scores for Groups 1 and 11 were 
30.5 and 29.0 which are much above the mean of the normally distributed 
scores of the achievement test. The maximum possible score for the achieve - 
raent test was 42. It indicates that both the groups learned matrix algebra 
very well. 

On the basis of their responses, the error rates were counted for 
each frame. Those frames which had higher error rates (more than 10%) or 
had ambiguity of language, were rewritten. 

Thus revised materials were given to other graduate students of 
the department of Education, twelve from the University of Chicago and 
seventeen from the University of Hawaii, to read, and the above procedure 
for revising the program was repeated. This time the students were asked 
to note down the total time taken to read the whole material* The average 
time to read the program is ten hours and thirty minutes. The revised 
program is given in the appendix. 
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INSTRUCTIONS 



Each numbered section in the text is called a "frame" and in most 



such frames a written response is called for at some point. Following 
the frame is the "answer" which is the correct answer. Therefore, 



proceed as follows: 

1. Use another sheet of paper to mark the printed "answer" while 
you read the frame and write your response in the blank space. 

2. Move the paper down and compare your response with the printed 
answer. If they are identical, proceed to the next frame. If they 
differ, mark your answer with a colored pencil mark, then proceed to 
the next frame. 



3. You can work according to your speed as long as you understand 



the material well 



Go to the next page and begin your Ttfork. 
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Unit I: Definition 

1. A matrix is basically a very simple way of organising the numbers used to 
describe the various kinds of information that come to our attention every 
day. For example the scores of two students in three subjects can be i7ritten 
conveniently as follows: 

Subjects 

Eng 

Students 



1st 


Eng 


Math 


Hist 


15 


50 


31 


2nd 

i 


14 


52 


35 



The score of the 1st student in History would be 31, and the score of 
the 2nd student in Hath would be 52* Thus a matrix is a rectangular array 
of elements. The above matrix has 2 rows and 3 columns and it is said to 
be a 2 x 3 (read 2 by 3) matrix. In the above arrangement, the numbers 15, 
50, 31 form the first roitf and the numbers ..., ..., ... form the second row. 
The numbers 15, 14 form the first column, numbers 50, 52 form the second 
column and the numbers ..., ... form the third column. 



14, 52, 35 
31, 35 
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6, If there are 4 numbers in each row and column of the array, the matrix is said 



to ^ ea by matrix or a matrix of size •««••••• 

(number) (number) (number) 



4 by 4 



4 



7. If there are n numbers in each row and column of the array, the matr 
said to be an by , , matrix or a matrix of size 



IX is 



n by n 
n 










8. The arrays 




x 





1 

| I 


1 




2 1 

i 


l i 

i . i 

1 and 


i 2 

( 


* i 

i 

* 


3 2 ! 


I 


1 


■ 3 


J 


1 ! 




j 


• matrices. 


, or matrices 


of 


size ... 



2x2 

2 



9. The array 



2 
- 1 



4 

0 



3 

1 



is a x matrix, because it has rows and 

(number) (number) (number) (number) 

columns . 



2x3 

2 (rows) , 3 (columns) 
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10 



f 



I 






t 






i 



t 
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10. The array 



3 

1 



4 

2 



is a ........ x 



0 -5 

matrix, because it has 



rows and ........ 



columns. 



3x2 
3, 2 



Comment 



Thus xtfe have seen that the order of a matrix is given by stating first the 
number of rows and then the number of columns in the matrix. If the number 



of rows is the same as the number of columns, then the matrix is square. 

Our idea is to consider such an array of many numbers as a single object . 

gjtray., a matrix and to give the whole array a single name or symbol. By 
regarding a rectangular array of numbers as constituting a single object, 
a matrix, we Ttfill be able to handle large sets of numbers as single units, 
thereby simplifying the statement of complicated relationships. Thus, we 



might call our matrices of frames 8, 9, and 10 as A, B, C, and D. 

1 “'~ MU, 



2 

3 



1 

2 



L 



2 

-1 



4 

0 



etc. 






I 



1 



i 



■5 



I 






I 



11# If E is the 4x5 matrix 



1 


2 


1 


9 


1 

10 ! 


! ° 


7 


8 


-2 


i 

12 j 


5 


-3 


4 


-5 


* 

15 


43 


17 


-10 


11 


22 



the element in the second rox -7 and third column of E is 8 and it can be 
written conveniently as [eJ ^ q 53 The element in the third row and 



r n 



fourth column is and it can be written as E j = -5, 

(number 



5 



3, 4 



12# In '.he matrix of frame 11, the element in the fourth row and fourth column 

is ••••••••••• 



11 



13# In the matrix of frame 11, [e],^ ^ “ 



» 11 



14. In the matrix of frame 11, | JS j 



4, 




li 



4 




I 



15. 



In the matrix of frame 11, the element of ith row and jth column can be 
written as [e | 



ra 






4 4 

J 

16, In the matrix of frame 11, there are entries in one rox7. 

(number) 

5 

17. In the matrix of frame 11, there are •»•••••• entries in one column. 





1 



13 



18, In the 4x5 matrix of frame 11, there are entries altogether. 



i 



20 



19. In a 3 x 3 matrix there are entries. 



D 



20. In a 3 x 4 matrix there are ••.•••••• entries. 



12 



21. In an n x m matrix there are entries. 



n m 












22. In an n x n matrix there are 



entries. 



a 



n 



n 



% 

i 
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23, Let I be the matrix 



0 



0 

0 

0 



0 

1 

0 

0 

0 



0 

0 

1 

0 

0 



0 

0 

0 

I 

0 



0 ; 
0 
0 
0 
1 



then 



(a) bl 4,4 = 






( b ) W 3,4 83 



n* 



0 



3 



t 



it 



1 






5 



1 








liitt— iiUM— — rtiimi —i nr — ‘i — -iir — • u— — rr- " - — 1 








^ /! --H 1 r^.-Vsv^ * 



R9HMI 
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(g) The element £l ] ^ j = when i ® j 



(h) The element [ij 






when i ^ j. 



(i) There are entries in the matrix. 



25 



l#*rf f 3 

£ 



24. 



If we interchange the rows and columns of a matrix then the new matrix formed 
will be different from the original matrix and will be called its transpose . 
If A is a matrix then A' will be its transpose. 

For example, if 





r* 




- 






* 


1 








i 


3 


2 






1 


2 






A = 








then A 1 = 


( 


3 


•*1 








2 


-1 


0 


} 




2 

» 


0 

M 






transpose of 
r 




















1 


3 


4 




• •• 


• 


III 


# 


MM 


B = 


2 


7 


2 


i.S B * 5=5 


« Ml 


1 1 I 


• 


• Ml 




0 


3 


-1 




III 


• 


III! 


1 • M 








MH 




Ml • 








- 
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Unit II: Equality of Matrices 



1, If 



r "i 


1 




4 3 1 


4 3 1 


2 5 2 


and B - 


2 5 2 


-7 3 0 




-7 3 0 









then 

(a) The size of A is to the size of B. 

(word) 



equal 



(b) All the entries of A are the 
of B. 



(word) 



as the corresponding entries 



same 



(c) The two matrices A and B are said to be equal if 
(i) A and B are of the same size. 

(ii) All the entries of A are the same as the corresponding 
entries of B. 

Therefore matrices A and B are 

(word) 



equal 



19 



2. Now if 



r 4 s ii 


i 


1 4 3 1 | 


i | 

A - | 2 5 2 

i.- 7 3 °. 


i i 

and B = 


i 

2 5 2 




-7 2 0 

• * .-w 



then 

(a) The size of A is «•••••••••• to the size of B. 



equal 



(b) The entry of row and column of matrix A is different from 

the corresponding entry of B. 



3, 2 



(c) Therefore all the entries of A are 
pending entries of B. 



as the corres- 

( same /not the same) 



not the same 



2J 



(d) Therefore the matrix A is 



to the matrix B. 

(equal/not equal) 



not equal 



3. If 



A 



4 

2 



5 

5 




then 

(a) The size of A is 



and B = 



to 



4 3 1 

2 5 2 

L 2 0 - 1 

the size of B. 



not equal 



(b) All the entries of A are •••••••• •••••••• •••••••« as the cor- 

responding entries of B, 



not the same 



(c) Therefore matrices A and B are 



not equal 



2i 
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4, Let A be an m x m matrix, let B be an n x n matrix, and suppose m ^ n. Then 




A is . to B. 

(equal/not equal) 



not equal 



5. Let A be an m x ra matrix, let B be an n x n matrix and suppose m - n, | aJ » 






i* |bJ for all i and j between 1 and n. Then A is to B. 

(equal/not equal) 



not equal 



6, Let A be an m x m matrix, let B be an n x n matrix, and suppose m - n, [aJ ^ j 

i* 

- [B J ^ j for all i and j between 1 and n. Then A is to B, 



(equal/not equal) 



equal 



7. Using the foregoing definition of equality, xne can express certain relation- 
ships more compactly, ----- For example, the equation. 



P- 




•* 


x + 2y a + 3b 




4 -3 


x - y a - b 




1 2 




j 


L J 



x + 2y = 4 

x - y = 1 

a + 3b * -3 

a - b =* ? 






4 



* 



i 



I 






?$ 






?? 



I 



1 





hi 
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S, Express the following relationships in matrix form 



X 


+ 


y 


— 


3 




y 


+ 


z 


= 


4 




z 


t 

r 


W 


= 


5 




*7 


+ 


X 


S5 


6 





1 



f 

j X + y 




i 






9. Let A = B 



where A 



then 



r 




** 


x + 2y a - 3b 


and B = 


0 1 


«x - y 2a + 4b 




% 35 






m* 




(a) x -1- 2y 









PS 







23 



(c) a - 3b = 



1 



(d) 2a + 4b = 



35 



10. Let 




0 0 0 
(a) [ B J ± 9 j = 0 when i , 



0 



1 

J 



(< or >) 



then 



J 



" ^ " (greater than) 



(b) 




1 when i 

(< orjC) 



m 

j 



ii 



it 



(less than or 
equal to) 



rv •* 

ZM* 






Unit III: Addition of Matrices 



| 

v 

I' 



1* To make the study of matrices meaningful and useful T\re must now consider 
basic operations with matrices. In this unit we define and study sums of 
matrices. Products will be considered later. 

Let the scores of Arithmetic and English of three students for the 
mid-term as well as for the final examination are represented by the 
following two matrices A and B, 

Arith Eng 





6 


8~ 




r 9 


n 1 
u 


A = 


8 


10 


B = 


9 


7 




7 


8^ 




8 


6 



Matrix A represents the scores of mid-term and matrix B represents the 
scores of the final examination. Now to find out the sum of scores on both 
the examinations for each student, matrices A and B will be added therefore 





6 


8 




9 


8 




6 + 9 


CO 

i 

co 


! 


fl5 


16 


A + B= 


8 


10 


d* 


9 


7 


=3 


8 d- 9 


10 d- 7 


l 


17 


17 




7 

mm 


8 

mm 




8 


6 




7 “!* 8 

MM* 


8 •!- 6 

+0 




15 


14 



The rule for the addition of matrices is: The matrices of the same size 

are added by adding each element of one matrix to the corresponding element 
of the other matrix. Thus, 








— J 



“ 

6 



2 



7 

5 



mmmmmmmatammmm 




' 










I 










\ 



■; 






3 







I 



1 





< ’ * " 

! 



PRPMiRilR 



5PSHWSB 



26 



;“ " i* 

i i 



t 

1 

• 1 


j G 9 




j 5 


16 


2. If A - | 


! 

11 13 

i 


• 

3 


B = ! 

12 


15 


of A and B, then: 




• 


- 



(a) 



C = 

*t=j 



and C represents the sura 



«» 



13 


25 


23 


28 



(b) The size of the matrix A is 



»•» 



2X2 or 2 



(c) The size of the matrix B is 



2 X 2 or 2 



(d) The size of the matrix C is 



2 X 2 or 2 



o 
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27 



3. Let 


**« 


1 ! 


r 1 


A = 


2 3 1 

-1 0 4 


1 

and B = 


-4 2 1 

1 3-2 




, «• 


i 


* 



and let C be the sum of A and B, so that 



C 11 


c 12 


°L3 






C 21 


c 22 


c 23 


* 

) 


then 



s 



(a) The value of c 



11 



v 



-2 



(b) The value of c 



22 



(c) The size of the matrices A and B is 



2 X 3 



I il KJC 










r 



n 









P* K. ♦ 



IER1C 



28 



(d) The size of the matrix C is 



2 X 3 



4 , The sum A + B of tv70 m x n matrices A and B is the m x n matrix C whose 

mm m 

entry in the ith row and jth column is the sum of j A J ^ j and JBJ ^ 



Therefore, 

(a) 



|a + B 






r -j 
A . . 









L B J 



i>j 



(b) t A ]i,j + Mi,. 



fA 



1 

+ Bl 



j, . . * 

• X »J (entry of matrix C) 



[cj • 4 

L 1 i »3 






i 



i 



i 



i 



aauca 




5. 






23 



If 



“ 

a, 

11 


a 12 


i 

i 


>u 


b 12 


j 

I 


C 11 


c 12 


a 21 


a 22 


X 

t 


b 21 


b 22 


= 


C 21 


c 22 


J*31 


a 32 




_ b Sl 


b 32_ 




C 31 


c 32 



then 

( a) c = 



a ll + b U 



(b) c = 



a 



12 



+ b 



12 



(c) 






21 






a 



21 



21 
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(b) Therefore A -5* 0 



(name of the matrix) 



A 



7. Let 





1 0 








i 

I 0 


0 


A = 


2 -1 


and 0 


ss 


0 


0 


t 


-3 4 

— 








0 


0 




*»* 

I 1 


o' 


i 


r° 


°i 




0 


= ! 2 
l 


-1 




0 


0 = 






[-3 


4 




0 


°J 


• f • • 

(wr 



then 



• • « • O 



(write entries) 



(name of the 
matrix) 




8. If A is an m x n matrix and 0 which is a zero matrix is also of size ra x n, 
then 



A + 0 - 




3 










32 




Comment 

TJe have seen that if a zero matrix is added to another matrix (of course, both 
of the same size) , then the resultant matrix is equal to the non-zero matrix. 
Therefore in matrix algebra zero matrix in addition is identity matrix as in 
the algebra of real numbers 0 is the identity element for addition (0 + a ** a) 




8 









h. 

\ l r- 






1 

4 



1 ^- 



3 

4 



and 



(b) B + A 



1 


i j 


F 




- 


I 

-2 1 


l 

t 


1 2 








+ 


| 






0 1 




4 3 












i 

1 








(C) 



Therefore C + D = 



D + C 



11. Let A and B be two matrices of size m x n so that 
A+B ■ C 
and B + A ~ D 
then 

(a) the entry in the ith row and jth column of (A 4- B) will be 





(b) and the entry in the ith row and jth column of (B + A) will be 



(c) The elements of a matrix are numbers, [a]. . and [b 1. . are elements 

' J i,j 

of ith row and jth column of matrices A and B, and therefore 
are 



numbers. 



(d) Numbers are commutative with respect to addition, as 2 *!■ 3 = 3 + 2 or 

more generally a + b = b + a. Therefore f A ] . . and IB 1. , are 

i Ji,j J x,j 



in addition. 



commutative 



(e) Therefore: 

and l B ii,j + i A ji,i are 



(equal or unequal) 



equal 






I 



I 









1 






i 



* 



I 



i 



4 










A 



o 

ERIC 






35 



(f) But 



Wi.j + Wi.j = H sj 



H,i + L A J*.j -L d Cj 

Therefore jc ^ and [d ^ are ..... 



and therefore C = D« 



equal 



(g) When C and D are equal then (A + B) and (B -!- A) are 



equal 



12, We have seen that matrices A + B - B + A and therefore matrices of the same 



size are commutative in addition. It is also true in the algebra 1 of real 
numbers where numbers are commutative in addition, e.g., a + b » 



b -!■ a 






aaat ^ aia 






mmm 
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Let r 


<■* 




mr 




1 




r 






' 1/2 


2 




-1/6 


4 




0 -1/3 




A = ■ 




B = 








and 


c = 








1 0 


-1 


9 




0 


-1/3 




1/4 1 




L 


J 




L. 




* 




L 


• 




then 




r* 




i 


r 


' 


! r 


- 






-1/6 i 




1 1 0 


-1/3 




f 

>•••••• • • 


» « t • • 


(a) B + C 


S3 






+ 


















O -1/3 




1/4 


1 




*••••• • < 


i • • • • 






+ * 






l 








- 



-1/6 11/3 

1/4 2/3 



and 


i 

t 


r 

t 


- 


*-* 

! i 


- 


! { 


i 




i 


1/2 


2 


-1/6 


11/3 




• • • 0 • • • • • • 


(b) 


A + B -5* C = A + (B + C) « 






X j 

1 | 


2/3 

•M 






0 


-1 


i 1 ' 4 




t • • • • •••«• 




1/3 17/3 

i 1/4 -1/3 




(o) Also r 

jl/2 

A + B = j 

1° 






1 

2 




-1/6 4 




+ 




1 




0 -1/3 

PM *■•«* 



I 

< 

cs t 

t 

I I • # • ^ 



1/3 

0 



6 







i 
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(d) And A + B -!- C = (A + B) + C = 



w ' J 






r 


j 1/3 6 

j j 


1 

+ 


0 


-1/3 




j 0 -4/3 j 




1/4 


1 


V 



• ••• • • • • 



• ••« • * a « 



1/3 

1/4 



17/3 



-1/3 



(e) Therefore 

A + (B J r C) and (A + B) + C are 



equal 



14. Let A, B, and C be m x n matrices. 
Then 

(a) f(A + B) + C | , = [A + Bl 

*. A *3 j > i 






+ fc) . * 



« • • • #••• •••• 



Wu'l'Lflj 



i ,J 






(b) And [A + (B +C)] = [A] + [ B cj . 

S J i J J 



+ 



Mi.J + f- C i , 



i»j ’* ~i,j 






•d 

I 



'l 



i 



V 
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Unit XV: Subtraction of Matrices 

1. Let A be an m x n matrix. Then the negative matrix of A is the matrix 
each of whose entries has the opposite sign of the corresponding entry 
of A. The negative matrix of A will be denoted by the symbol -A. That 
is, the negative matrix -A is defined by - : A . . = *-A| . .. 

If A and B are matrices of the same size, then the difference of A and B, 
denoted by A-B, is the sum of A and the negative of B. That is, A-B is 
defined as A -!■ (-B). 

If 

A = -1/2 A 

0 1 



then 











t 

11 

<J 

I 

/N 

CO 

vy 


<3- 

CM 

rH 

1 








0 1 















1/2 -4 

0 — 1 



-1/2 4 


1 

» 


1/2 -4" 






o 

J 




t 

O 

/ 




MMI ••••• 

•M «*"> 



(b) 

A *!* (-A) = 







■^rtfsss«r*t^«X3c^ 
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i 





(b) 




r 




i I 


(N 

t 








-1/2 


4 


! 


2 i 


(“A) + (-B) = < 


; - 






. j- * 
/ ' > 

j 


t " i 




| 


| 


0 


1 


1 1* 


2 i 

J 





• • • • • • • • 




• • • . •••• 




**" • • • • • • • %* 






+ 




rr 






• * • • • • • • 




•••• •••• 


1 


j • • • • • • • j 



1/2 


-4 


f -2 
* { 


- 2 l 


; 1 
J 

j 


- 3/2 


- 6 * 


0 


-1 , 


u 


-2 1 
J 


= 1 

4 

: 1 


1.-1 


-3 



(c) From the results of (a) and (b) what relation can be established between 
-(A 4- B) and (-A) * (-B)? 



-(A B) = (-A) + (-B) 



4* If A is an m x n matrix, then 






A , -A ; A t A ; 0 



5, In the above frames we have seen that 



A (-A) = 0 



It is true for any size of A. 



43 



6. If 



x 3 


2y - 


wm 

s 




0 


-6” 


a 1 


4x 


6 


= 


-3 


2x 


b-3 


3b 






2b 4 

IT- 


-21 



then 

(a) The value of x that satisfies the matrix relationship is 



-3 



(b) The value of y is 



1 




(c) The value of a is 



-4 



(d) The value of b is 



m 






* / v? s r ? t -f ‘fHff-y 
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8, Compute 



1 1/2 

i 


1/3; 

i i 


1/6 


1/7 


1 

f 1/4 


i' 

1/5 5 


j 1/8 


1/9 



2/3 


10/21 


3/8 


14/45 



i 



9 • Does the sum 

I 3 

! 

i i 
j 3 

make sense? 



2 

3 

1 



0 

0 



o i 



o 



No 



10, Does the sura 

: 3 2 i ; 

t 

t 

j » 

; 1 3 2 ! 

i i 

1 3 1 2 : 

! j 

make sense? 



+ 



; o 

\ 

; o 

i 

; o 



o 

o 

o 



0 : 



0 i 



0 




Yes 






i 



I 



I 

i 






mm m 
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11, Compute 



1 


2 


o \ 

° s 

i 


9 

\ 


8 


7 


i o 
» 


0 


o ; 

1 


. 10 
1 


10 


10 


4 


5 


6 i J r 
< 

| 


j 6 


5 


4 ' - 

r 

t 


i 

1 0 

i 

l 


0 


o ! - 


! 10 
i 


10 


10 


7 


8 


9 ! 


1 3 


2 


1 : 


j 

: 0 


0 


* 

1 ! 


• 10 


10 


10 



i 



to 0 o 1 

j 1 

jo 0 0 j 

; I 

; o o -l ; 



(Comment) 

Insofar as only addition and subtraction are involved, the algebra of 
matrices is exactly like the ordinary algebra of numbers. 

Suppose A and B are know* matrices of the same size. Now consider the 
matrix equation X A ® B, where X is of the same size as A and B but is 
unknown, If this equation involved ordinary numbers rather than matrices 
we would simply add -A to both sides of the equation; i.e. , 

X -! A (-A) = B (~A) 
or X = B - A 

It happens that x$e can do the very same sort of thing with matrices. 



o 
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12. If C and D are knom matrices of the same size and Y < C = D, where Y is 
of the same size as C and D but unknown, 
then Y = - 






for the matrix X. 




o 












1,111 W.l III C ,» . ,1 II n I IPPP! R . I III ,IUI 1 1 1 l.,l| Jl |l. H* WW1"WI I III 1 ; 1 









i 




49 



Unit V: Numerical Multiples of Matrices 



1, Once we know how to add numbers, it is customary to define 2x as the sum 
x -i- x, 3x as the sum 2x x, etc. Fractional parts of x are defined by 
requiring that (l/2)x+ (l/2)x = x, (l/3)x (l/3)x (l/3)x - x, etc. All 

of this can readily be done with matrices. If we add two equal matrices, 
the sum is clearly a matrix in which each entry is exactly twice the 
corresponding entry in the two given matrices. Thus 



! 2 



i 



|-1 

1 L 



■M 

\ 

. | 


** 2 


3 


4 


4 


6" 




2(2) 


2(3) 


) 


-1 

•* 


0 




-2 


0 




2(-l) 

■M 


2(0) 



The sum 3A = 2A -!• A = A + A *!• A is equally clearly the matrix each of whose 
entries is exactly three times the corresponding entry in A. The equation 
(1/2)A + (1/2)A = A defining the matrix (1/2)A is clearly satisfied by 

the matrix each of whose entries is exactly one-half the corresponding 
entry of A. 



2 , Compute 



2 


•* 

3 


r 


-1 


0 

** 





6 9 

-3 0 



1 

- 1/2 



3/2 

0 












3 



* 



& 






1 



i 



i 



3 



% 



Compute 


2 3 


1 


■* * 




1/2 




= 








-1 0 










L *5 




JWf ^ 4 



1 



1 

i 























5'j 



4. 



Let A be an n by n matrix and x any real number. Then the product of A by 
the number x is the matrix whose entry in the ith row and the jth column 
is x times the corresponding entry {*] ^ j in the matrix A, That is, the 



matrix xA is defined by the formula 

Hu - * 



Hu 



Notice that the product of a matrix by a number is another matrix. 



5, Let x = 5, y = 2 and A = 
then 



1 3 

2 1 



(a) x(yA) « 5 S 2 




’] ■ I.' ] 



2 


6~ 




10 


30 


4 


2_ 




20 

m 


10 



(b) and 

(xy)A =(5x2) 



1 3 

2 1 



= 10 



1 

2 



3 

1 



a 



10 


30 


20 


10 
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i 



t 
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(c) From the results of Frames 5(a) and 5(b), x(yA) and (xy)A are 




equal 



6. If A be an m x n matrix the relation x(yA) - (xy)A will still be true as 
is clear from the following 



[ 



x(yA) 



i 






x yA . . 

13L»J 



= X 



<yj A j 



i 



xy [ A Jij = Ny 



7. Again let x = 5, y = 2 and A = 



1 

2 



3 

1 



then 



(a) 


'1 3 “ 




1 3 




• • • • 


(x y)A = (5 + 2) 




* 7 










2 1 




to 

1 




• ••• • t • • 



7 21 

14 7 














Si 



52 



(b) And 

xA -!* yA = 5 





t 






r 7 


21 ' 


*r 

J 




r* 


j 14 


7 

. ^ 



5 


«• 

15 




“2 


6 " 


10 


5 _ 


i 


4 


9 

> 



(c) From the results of Frames 7(a) and 7(b) we have 

(3C *: y)A = (xA) •: (yA) 

and it is true also when A is an m x n matrix. 



i c. if 

* 


l 

A * 


* 

\ 

1 


3 


i 

s 




2 


1 



then 



1 


1 3 






(a) (-l)A = (-1) 




r - 




i 


2 1 




W 



! 

I. 



I; 



l 

I 

f * 



~ 1 -3 

-2 -1 






\ 



! 




mmimimmmimmmimmimmimmmmiimmmmimmm 







mmmm 
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and 

(b) 



xA + xB.= 5 





1 : 


r 

1 


1 


! 1 3 


l 

1 v 5 


1 

O 


4 


2 1 


i • , | 


-1 


0! 


j 










" 




f m 1 


r I 


1 1 • • • • • • 


i 

j ] 

1 


j#### • •#«.« j 

= 1 


[»••• 

j 


• • • • • « ♦ • 




i j 

• ••« • • • • j 


• • • • 


V- J 


1 1 


L .? ] 


L -J 



f 1 




-1 


l 


j 5 15 




10 20 I 




» 


^ 1 

*» 




1 — 


1 10 5 

1 




I 

Ui 

0 




j 






L 



15 

5 
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(c) Therefore from 9(a) and 9(b) we have 
x(Av-B) = xA + xB 

and it is also true when A and B are m x n matrices as may be seen from 
the following: 

[■ 

= X S 



x(A B)j 

* 1 jj 



1 



X, 



= xi A Bj . . = xsj A ! . . *•'{ B j. . > 

l Ji,j y- L ,J j 

= x( A ‘ x[b i 

- -'i.J 1 'i,j 

r 



Mi,] + i xB Ji,j 

[kA •!- XSj, 



10. Also if x = 5, then ( 



x • O 2 = 5 



r- 

0 


0 

1 






-1 


0 


0 









0, 



0 

0 



0 

0 




11. It is true in general that any numerical multiple of the zero matrix is 
the zero matrix, whatever the size of the zero matrix may be. It can be 
proved easily by using our notation for the entries in a matrix as follows 

[ x • °"3 i,j = x C°*3 i.j = * • 0 = 0 



Therefore x * 0 =0 

n n 



12. When 




then 



0 • 




3 

1 




HHfc 

0 0 

0 0 



13. The product of any matrix by the number zero is the zero matrix, whatever 
the size of the matrix may be. It is seen from the following: 

Using our notation for the entries in a matrix, we may simply write: 




Thus 0 • A = 0 n 



i 0 

; ERIC 









iipp» 



WSSsi 



14. Let 



A = 



then 

(a) 



3A = 



and 2C - - 



55 



) 2 1 -3 ' 


i 

i ! 


’ 3 0 5 | 




5 -1 O’ 


f 

jl 0 4 

*• i 


, B = 

1 


6 9 -1 j 


, and C - 

1 1 


7 3-1 

1. 



4B = 




(b) 





6 


3 


-9 ! 


i 12 


0 


20 j 


r io 


3A = 








1 ; 4B = ! 

! 24 






2C = j 




3 


0 "• 


12 


36 


1 14 








1 


1 L 




J 


L 



4b 2C = 



© • • • » 



-2 

16 



0 



-2 



22 


-2 


20 j 

I 


38 


52 


I 

-6 | 



| 










58 
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Unit VI: Multiplication of Matrices 



So far we have defined and studied: 

(1) Addition of matrices 

(2) Subtraction of matrices 

(3) Multiplication of a matrix by a number 

How we shall define and study the product of two matrices. 

The definition of multiplication of two matrices will be clear from the 
following example: 

1. Suppose in a University three departments— Education, History, and 
English— require teaching faculty in three categories— Professors , 
Associate Professors, and Assistant Professors. The following is 
the number of such personnel required in each department: 



Departments 

Education 



Associate Assistant 
Professors Professors Professors 



History 

English 



10 

5 

8 



12 

8 

15 



15 

10 

20 



We can, if we like, arrange this table in the form of a matrix: 



Personnel 



Departments 



| Education 

4 

{ History 
English 



t. 



CO 


03 


03 


u 


<U U 


4J M 


o 


•U O 


CJ O 


CO 


« 03 


Cfl 03 


CO 


•i4 03 


4J 0) 


0) 


O 0J 


03 <1) 


M-i 


O *4-4 


•H 44 


o 


03 O 


03 O 


u 


CQ 54 


03 54 


PW 


< PM 


<1 CM 


10 


12 


15 


5 


8 


10 


8 


15 


20 



= A 















1 















I 



r. 






\ 



Suppose now the College has to pay the following money (in thousand dollars 

■i 

for each personnel: 



Regular Salary Summer Salary Transportation 



professor 




20 


4 


3 


Associate 


Professor 


15 


3 


2 


Assistant 


Professor 


12 


2 


1 


Again, we can 


arrange our 


informat ion 


in the form of a matrix: 





Expenditure 



G '• 
O 

4-J 

cU 



U 





Regular 

Salary 


Summer 

Salary 


o 

Cu 

C/3 

a 

u 

EH 


Professor 


~ 20 


4 


3 ~ 


v Assoc, Professor 


15 


3 


2 


Assist. Professor 


12 


2 


1 



= B 



It is clear that by combining the above two sets of informat ion, we can easily 
compute the expenditure of each department for the three categories of expendi- 
ture. That is, by correctly combining the entries in the two matrices, we 
can figure out another matrix which will have only departments and expenditures, 
(a) For instance, to compute the expenditure of regular salaries for the 
Department of Education, we figure as follows: 

Department of Education will spend for the regular salaries for: 
10 Professors at the rate of 20 thousand dollars for one + 12 
Associate Professors at the rate of 15 thousand dollars for one 
+ 15 Assistant Professors at the rate of 12 thousand dollars for one 



= 10 x 20 + 12 x 15 + 15 x 12 = 560 thousand dollars 



j o 
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We multiply each entry in the Education row of the matrix A (reading from 
left to right) by the corresponding element in the Regular Salary column of 
the matrix B (reading from top to bottom) and then added all the resulting 

products. This procedure is perfectly general. 

(b) Thus, to find the expenditure of summer salaries for the Depart.ment 
of Education, we multiply each element in the Educatio n row of the 
matrix A by the corresponding element in the Summer Salary column 
of the matrix B and then add, getting a total of 

10 x 4 12 x 3 15 x 2 = 106 thousand dollars 



Look at the example above and ansvjer the following: 

(c) To find the expenditure of the Department of Education fo l paying 
transportation, multiply each element in the Educat ion row of the 
matrix A by the corresponding element in the Transportatio n column 



of the matrix B and then add, getting a total of 

v x *.=69 thousand dollars 



10 x 3 12 x 2 + 15 x 1 



(d) To find the expenditure of the History Department for regular salaries, 

we multiply each element in the row of the matrix A by the 

corresponding element in the -------column of the 

matrix B and then add, getting a total of 

x x -b =340 thousand dollars 



History; regular salary 
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(e) Find the expenditure of the History Department for summer salaries. 



x 0 x x =,... .thousand dollars 



5 x 4 + 8 x 3 •!- 10 x 2 = 64 



(f) Find the expenditure of the History Department for transportation. 

-1- x. . . . ,-i- x. . . . ,= thousand dollars 



5x3 3x2 10 x 1 = 41 



(g) The values xdxich vie have in (a) to (f) are some of the elements of 
the matrix x-jhich is the product of matrices L and B. Ue can x-jrite 
these matrices as folloX'js: 



•f 

I 



ilO 

1 

j 


12 


i 

15 


t 


I 20 4 

i 


3 j 


5 


3 


10 


X 

i 


1 15 3 

j 


2 


i 8 


15 


20 


l 

f . 
1 
l 


j 12 2 


i 

I \ 


L 




i 


1 


L 


C 

* 


10 


x 4 + 


12 x 3 *! 


■ 13 x 2 = 


106 


5 


x 4 -!- 


8 x 


3 


10 x 2 = 


64 



k 

I 



10 x 20-1- 12 x 15 + 15 x 12 = 560 = a 
= 5 x 20 -!- 8 x 15 -!* 10 x 12 = 340 = d 

L h 

b 10 x 3 12 x 2 -1- 15 x 1 = 69 = c 

e 5 x 3 *!- 3 x 2 10 x 1 = 41 = f 

i j 



I 

*5 



| 










mmrnmmm 



tMsm 



2. In the above frame the value of h can be obtained by multiplying each element 

i n the row of the first matrix by the corresponding element in the 

(which row) 

column of the second matrix and then by adding. 



(which column) 







third 

first 




3. In the matrix above [frame 


1(8)], 






(a) the value of i is 








. . . , .x *iv. , . 




..x =117 








8 x 4 15 x 3 *!- 20 x 2 




(b) the value of j is 


. , . . .X v. . , . 









0x3 15 x 2 20 x 1 = 74 



COMMENT 

The most concise description of the process of multiplication is: Mu It ip ly 

row by column/ 1 Very simply the rule is to multiply entries of a row by 



corresponding entries of a column and then add the products, Thus., given 



I 






f 



I 



t 



I 



two matrices A and B, to find the entry in the ith row and jth column of the 
product matri x AB, multiply each entry in the ith row of the left-hand factor 
A by the corresponding entry in the jth column of the right-hand factor B, and 
then add all the resulting terms. If A and B are 4 x 4 matrices, then the 
entry in the ith row and jth column of the combination matrix AB will be 

I M i,l ! B 1l,j i A | i,2 [ B J 2,J • ; [ A J 1,3 L B b,j i A ii,4 i B h.j 



4. 



5. 



A 



B 



AB 



1 2 




5 0 | 

! = 

I 


3 4 

J 


i : 


6 7 ! 

! ! 



(3 x 



5) + 


(2 


x 6) 


a 


X 0 + 2x7) 


5) + 


(4 


x 6) 


( 


,...x ) + (.... 




write numbers 



,x. 






0 + 4x7 



A 


B 




1 

o 

r-A 




r 4 3 2* 1 


| (4 + 0 -!- 

i 


2 3 4 




7 0-6 




(8 + 21 • 


-4 5 -2 




-3-5 6 j 


j (-16 +35 

r x 


J i 


«•» ^ ^ 





AB 



(write numbers) 



30 - 12). 



(write numbers) 



2 + 0-6 



-12 o + 10 
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I 

% 

| 



1 



1 

I 

8 

1 



1 

! 

1 

% 

1 

I 



1 



« 






I 



I 



I 



1 



8 



I 



:?■ 

1 



! 



.i 
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6. Let us try the following schematic device to explain multiplication : 



B 



5 



•2 ' 



* / I 

. / Ji i ) 

V/ / *• V - 

y ''/ •dVx' 

% . ' y y / 

/ <o / , 

/ 



/ V +/ 
/ \ / 



a / /' y 

/ *' * 

6 ) l - J : 

3 4 



6 i 

■i 

9 ! 



AB 



I r; 



p x 1 6 x 2 *!~ 7x1 



2 x 1 + 3 x 2 4x1 



5 x 6 6 x 5 -I- 7 x 9 

. , , x. . , ,X. . . ,x. 

(write numbers) 



2 x & 3 x 5 -!- 4x9 

We see that the matrix A is of size 2x3 while the size of B is 3 x 2. Though 
the matrices A and B are not equal, yet there is an entry in each rox-7 of the 
matrix A to match with each entry in a column of the matrix B, and conversely. 
It follows that the product is not defined unless the number of columns in 
the A matrix is equal to the number of rows in the B matrix. When the number 
of columns in the left-hand factor equals the number of rows in the right-hand 
factor, the matrices are conformable for multiplication . 
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8. Perform the follov7ing matrix multiplication: 



1 2 3 j 


1 S 1 0 1 ! 


1 


j 


4 5 6 | 


2 1 2 


11 


. - 




| 






7 8 9 


! 4 1 3 




i 

| 


•r* -* 


' “ - 




.i 



17 


5 


14 ! 


38 


11 


32 


59 


17 


50 



9. Perform the follov7ing matrix multiplication: 




3 1 

-1 2 



1 

4 



2 

■* 

0 



3 

1 




iU. Perform the following matrix multiplication: 



1 


0 


o" 


1 

j 


~ X 1 


x 2 


x 3 






i • o • • • « 


***1 


0 


i 


0 




y l 


y 2 


y 3 


= 








i- 1 

O 


0 


1 


1 


Jl 


z 2 


z 3 




* 4 







in « 




















Unit VII: Properties of Matrix Multiplication 



1 » Let 





0 


"n 

o j 


0 


i 

1 


V = 




! and B - 






i 

i 

I 




0 


0 


0 


1 




? 







then , 



(a) 

AB = 


0 


0 | 


l 

1 


0 


1 

JL 


- 






j 

1 




1 


0 




0 


0 

1 


t 



w 




■o 



tf 

I 

S 

I 

$>■ 

| 

| 

I 

I 



I 



’ 0 0 

0 1 j 



and 

(b) BA = 



I 0 1 




0 0 




( 


[ 0 0 




1 0 

1 . J 




1 

i 



a 



, 



i 



1 



Oi 



0 Of 



In the previous frame we saw that 

AB and BA are 

(equal/not equal) 



jo 


0 




1 


■n 

0 


AB = 




and BA = j 


l 

1 




0 


1 


1 

i 


i 0 


0 



therefore 



(not equal) 



, 0 
ERIC 










} 



then 




s 



• c 



\ 




i 

v' 







9 2 5 

7 6 10 

7 -2 -1 

t 



i 



and 

(b) BA = 



I 



4 



2 3~ 

0 1 






o 

ERIC 



{yZBBfflfl333 



i 4 
3 



10 " 

10 



6. Again we saw that 



AB = 



i 9 


2 


5 














4 


10 


7 


6 


10 


, and BA - 














3 


10 


7 

•* 


-2 


-1 




l - 


.J 



, therefore AB and BA 



are .matrices 

(equal /different) 


















i 



I 



I 









i 

i 

i 
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I 

| 



different (or not equal) 



COMMENT 

We have seen that if there are two matrices A and B which are conformable for 
multiplication, then their product AB is not equal to BA. Thus we have a 
first difference between matrix algebra and ordinary algebra, and a very 
significant difference it is indeed. When we multiply real numbers, we 
can rearrange factors since the commutative law holds, e.g. , for all values 
of x and y, xy = yx. When multiplying matrices, we have no such law. We 
must, consequently distinguish between the result of multiplying B on the 
right by A to get BA and the result of multiplying B on the left by A to 
get AB. In the algebra of numbers, these two operations of "right multi- 
plication” and ’'left multiplication” are the same; in matrix algebra, they 
are entirely different. Another way of putting this is: in general 

multiplication of two matrices is not commutative. 



I 

I 



I 

X 

I 

i 

1 






i 



i 



j 

| 



>* 

I 

I 

! 

I 

X 

I 

; 



I 
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COMMENT 

Vie have seen that the product of two matrices can be zero without either 
of the two matrices being zero. This is a second major difference between 
ordinary algebra and matrix algebra. In ordinary algebra xy = 0 only if 
either x or y is zero, but it is not true in matrix algebra because AB can 
be zero even if neither A nor B is zero. 



10. Now let 



A 



1 


2 


0 


( 


1 


2 


3 


and | 

i 1 


i i 

I 


2 


3 


1 


1 


0 


; • b = 


1 


1 


-1 


c = 


1 


1 


-1 


-1 


4 


0 




2 


2 


2 


i 


1 

L 


1 


1 



then 



(a) 


' 1 2 0" 




[ 1 2 3 | 

1 

t 


AB = 


110 




1 1-1 




-1 4 0 




2 2 2 

. . 1 



4 



*3 



3 2 

2 ~7 



and \ 


| 1 2 0 | 


\ \ 


I 1 2 3 




1 


(b) AC = 


1 1 0 | 


1 ' 
P 


11-1 








-1 4 0 




1 1 1 j 


1 1 


— - 




5 

\ 



1 



A 

■5 



l 



| 



1 






I 



i 



> 

j 



£ 



1 

3 



I 

I 

% 






1 



3 



j 



l 



i 



.1 



I 



I 
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s 

jjj 


3 


4 


b - 1 


V 

I 

h 


2 


3 


2 


1 


3 


2 


1 

1 



\ 






l 



t ~ ' 



i T*» 



% 






o 
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11. In the previous frame we found that 





CO 

b-» 

J 




3 4 1 


AB = 


2 3 2 


, and AC - 


2 3 2 


4 


3 2-7 




3 2-7 



therefore AB and AC are 



(equal/not equal) 



equal 



Here we sec; that AB - AC, while A f- 0, and B /- C , while in ordinary algebra 
if ab = ac then either a = 0 or b - c (cancellation law for multiplication). 



Therefore, the breakdown for matrix algebra of the law that xy = yx and of 
the law that xy = 0 only if either x or y is zero causes additional differences. 
The cancellation law for multiplication for ordinary algebra can be proved 
as follows: 

(a) ab - ac 

(b) ab - ac = 0, 

(c) a(b-c) = 0, 

(d) b-c = 0, 

(e) b = c. 



\ c / u • 

For matrices, the above step from (c) to (d) fails and proof is not valid 



for matrices. 



iiBiiaiHBaiiaaaiiaiiiiBfiiliaiiM 





COM M E H T 



I 



f 



I: 



Let us consider another difference. We know that a real number a can have 
at most two square roots; that is, there are at most two roots of the 

* 

equation xx = a. 

Proof, Again, we give the simple steps of the proof: 



(a) 


Suppose that yy - 


a ; then 


(b) 


xx = yy, 






(c) 


xx- yy = 0, 






(d) 


(x-y) (x-l-y) = xx -! 


• (-yx 


+ xy) - yy. 


(e) 


• 

ii 






(f) 


From (d) and (e), 


(x-y) 


(x-: y) = xx- yy. 


(g) 


From (c) and (f ) , 


(x-y) 


(x-’.-y) = 0. , 


(b) 


Therefore, either 


u 

i* 

>> 

i 


0 or x -!- y = 0 


(i) 


Therefore, either 


x = y 


or x = -y. 



For matrices in general, statement (e) is false, and therefore the steps to 
(f) and (g) are invalid. Even if (g) were valid, the step from (g) to (h) 
fails. So the proof is completely wrong if we try to apply it to matrices. 
In fact, it is false that a matrix can have at most two square roots. 



12. Multiply the following: 



X 

o 

i 




r* 

0 x 






1/x 0 

I ; 




1/x 0 


i 


• ••<*• • • • a • 

L 



1 

0 



o| 



o 

iERIC 
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Therefore , 



0 

1/x 



x 



0 



and 



0 x 

1/x 0 



are the square roots of 



1 

0 



If we give different values (except 0) to x, we can get different square 
roots of ! . Thus the verv simDle 2x2 matrix has 



1 0 


1 ' 


1 0 




. Thus the very simple 2x2 matrix 




0 1 




0 X 



infinitely many distinct square roots. 



13. Let A and B be two matrices of equal size, then 

2 2,2 2 

(A -1- B) (A J .~ B) = A + AB BA -!- B /= A -I- 2AB + B , 



because in general AB and BA are 



(equal/not equal) 



not equal 



14. Let A and B be two matrices of equal size, then 

2 2 

(A J r B) (A - B) = A - AB v BA - B ^ A^ - B^ , because in general AB 
and BA are 



(equal/not equal) 



not equal 



15. 



Let 

A 



1 


-1 


and B = 


1 


0 


0 


2 




1 


2 


„ 











then 



(a) the sum of A and B is 





1 


-1 




1 


0 






A + B = 






+ 




= 






0 


2 




1 


2 





















o 

ERIC 
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79 



I 

■H 



l 







(h) Now 



(A 



j. 



therefore 



! 


r 3 -6 1 


i i 

i 2 2 


2 -7 


B) (A B) = 




| , and A -I- 2AB -!- B = 


r 




6 15 | 

L ~ 


1 

1 


7 16 



2 2 

(A + B) (A -!• B) is to A 2AB -!- B 



(equal /not equal) 



not equal j 

! 
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* 

not equal 



16. Let 

I 0 0 

A - ] 

1 0 



then 



2 

A 



0 


o ! 


! So 


°i - 




1 


0 

i 

i 


i 

! * 


o! 

' 


► 



0 

0 




0 



17. Let 



A = 



0 

0 



n 

0 

. 



(n 



then 



A 



2 



j, 0 n 

-0 0 . 



'■J 



0 



any real number) 



oj 



fo 

0 



0 

0 



0 



18. Let 



j 

i 


! 1 2 




i o] 


A = 


I 

1 


and B = 




! 


3 4 




-1 1 



then 


1 


z'i 


" i o] 






(a) AB = 


3 

•* 


J 


-1 lj 














mmm 







I 



1 



ti- 



nt 









EJ 



f 

I 



I 



I 



L 



I 



S’ 



O 
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19. Let 

A 



I 1 



3 



2; 



f 1 



= ' 4 



"I 

1 



6 •: 



and I = 



0 



t 

* 5 



1 



2 ; 

.i 



! o 



0 

1 
0 



0 



0 



then 

(a) 



! 1 



AI = 14 



: 5 



1 

1 



2 

6 

2 



i 1 



! 



5 0 

i 

i o 



0 

1 

o 



0 

0 

1 



= A 



i 1 

! 4 
5 



1 6j 



2 | 



and 

(b) IA = 



! 1 


0 


oi 


! 1 


3 


CM 




o 


1 


j 

oj 


1 

t 

j 

I 4 


1 


6 




i 


! o 


0 


i < 

l 

1.- 


i 

t 

i 

l 

i 5 


1 


aj 


i 

i 





(name of 
the 

matrix) 



; l 

4 



! 



! 5 



3 

1 

1 



2 

6 

2 



A 



’I 



1 



1 



i 



! 



£ 

i] 



-i 



| 

:s 



1 



J 






I 



5 






i 






Therefore 



fil 



% 



i 



i. 



f. 



r. 



| 



& 



i 












Pr 

] 



o o 
oj 



(c) 



IA = A = AI 



20. (a) Let 



A = 



1 

2 

3 



4 } 

t 

5 i and I 



1 

0 



0 

i! 



then 



AI * 



i 1 

i 

l 

2 

1.3 



1 

f 


u 


0! j 


i 


i 


* I 

M ! 




0 


1 




j . 1 



< • ••<"• 



(name of the matrix) 



1 

2 



3 



4 | 

5 ! 

\ 

1 

6 i 



A 



(b) Now let 



I = 



i 1 
1 


0 


01 

i 


1 1 
j 1 


4 


0 


I 


o 

1 


! and A = ! 2 

' i 


5 


0 


0 


i 

i 

i 


I 3 


6! 



then 



1 
2 
1 3 



4 i 

5 = A 



6 1 



ERIC 
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i 



s 






is 



% 



1 



3 



4 



s 



1 


il 


0 


0 

i 

j 


t < 

j 


i 1 


4 




• •• 0 # * * * * # t 

L 


I 

IA = 


1 

i 0 


1 


0 


i t 

1 1 

1 


2 


5 




J 

fc 


3 

i 

1 


t 




1 

1 


| 

( 






i 


\ (name of the matrix) 


i 

i 


0 


o 


X 

..1 


i 


3 


6 ! 

. 3 


[ 

1 


s 



I 



I 



'■ 



I 



J 






rnmmmmmmmmmmk 










(c) Therefore 






r 



jt'f 

% 



4 



r 



I 



S 



§ 



ri ■ 
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AI = A - IA 



COMMENT 



We have seen above that the matrix I when multiplied to the matrix A (either 
on the left or on the right) just gives A back again. The matrix I is called 
the unit matrix, or the identity matrix for multiplication. The unit matrix 
for the multiplication of matrices plays the same role as the number 1 does 
in the multiplication of real numbers. (For all real numbers a, la — a = al.) 



The matrix 



*2 ■=! 



! l 

o 



0 



1 I 



is called the unit matrix of size 2. The matrix 



h ■ 



0 0 
1 0 
0 1 



is called the unit matrix of size 3, and so forth. 

i 

In general, the unit matrix of size n, which we may denote by the symbol I 
is the matrix whose entries are 



I . . = 0 if i j 

ns i,j J 



j I n j i,j = 1 if i = j , that is j I j = 1 for every j. 

Now let A be an n x n matrix. Then by definition of multiplication, the 
entry in the ith row and jth column of the product AI is 

L A ]i ,1 C 1 1 1 , j ' : i A Ji,2 M 2 ,j ' ! l A li,n 

Since , 1 ik, j - 0 whenever k is different from j , every term but one in this 
last expression is equal to zero and drops out, and we are left with just 



■ n »J 



1 



Thus 



one term: [ Aj^. fl] j>;j , since | f| ^ -l.f Aj^ fljjj 

the entry in the ith row and jth column of the product AI is simply { A I . 

‘ .1 » J 

Thus AI and A have exactly the same entries, and therefore AI = A. 
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Similarly 

Hu 



Wi,l! A a,J - L 1 ! i,2 : A ' 2) - 

i 1 1 . . : A:. . ; Since 5 . , =0 when k / i 

5. .11,1 . tl,j » i.k 



4- 1 T ! 1 A : 

••••••••«• * jx # . ,n. • 

’■ n,j 



= ! A i . . ; Since 1 1 ! • • = 



— 1 



t 1 1,3 



Thus I A = A 






and therefore 


A I = 


A = 


IA. 




Compute 














0 


0 


"f 


0 


0 


n 




0 


1 


0 \ 

I 


0 


1 


0 




l/n 


0 


1 ! 

0 \ ! 

t > 


; l/n 


0 


0 



t 

l 



(name of the 
matrix) 



I 1 

1 0 

I 

i.o 



0 

1 
0 



0 ; 



0 ; = 1 



22. In the previous frame we saw that 



i 0 
0 



[ l/n 
therefo re 



0 

1 

0 



n 
0 

0 : 
0 

0 



0 

0 

l/n 

0 



S l/n 0 



n 



0 f 



0 



0 n 

1 0 1= I 

l 

} 

0 0 1 



is the 



of I 



square root 



Notice that giving infinitely many values (except 0) to n, I can have 
infinitely many square roots. 
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23. Let 


•* 

4 


5 


6 ] 




1 


0 


0 


A = 


1 


2 


3 | 


and I = 


0 


1 


0 




L 7 


8 


9 j 




0 


0 


1 



then 



(a) 


4 5 6! 




1 0 0 




4 f t ft#tt -rtrr* 


AI = 


1 2 3 




0 10 


:= 






7 8 9 


, 1 


0 0 1 







(b) 

(AI)I = 



t 



(c) 

I • 3>-f 0 

f 

0 



(name of 
the matrix) 



4 

1 

7 



5 . 6 

2 : 



8 



A 



f 4 


5 


6] 


t 

f 


1 


0 


O 






1 


2 


3 


l 

i 


0 


1 


0 


=s 




1.7 


8 


9 I 


i 

| 

t 


.0 


0 


1 







(name of 
the matrix) 



4 

1 

7 



5 

2 

8 



6 

3 

9 



= A 



0 

X 

0 



0 

0 

1 



»♦» 


1 


0 


0 I 








0 


1 


0 








0 


0 


1 







(name 
of the 
matrix) 



l 



I 



i 



'i 



% 



2 



y 









I 






i 



I 






n 



I 



5 





f 



(e) 


4 5 6 




1 0 0 






((AI)I)I = 


12 3 




0 10 


■ := 






7 C 9 




0 0 1 







l 

'( 

I 



(name of 
the matrix) 



f 

% 



j? 

% 

fi- 

fe 

i 

t 

h 

! 

% 



4 

1 



5 

2 



6 

3 



A 
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Unit VIII: The Laws of. Matrix Multiplication 

We saw in the previous unit that two basic laws which govern multiplication 
in the algebra of ordinary numbers break down when it comes to matrices. 
That is, we are faced with the 

Breakdown of the commutative law : The product AB of two matrices may be 

entirely different from the product BA of the same two matrices. 

Breakdown of the law of cancellation : The product AB of two matrices may 

be zero even if both factors A and B are other than zero. 

Aside from these two laws, most of the other basic laws of ordinary algebra 
remain valid for matrices. These laws will be stated in this unit. 



1 , 



Let 

A = 

then 

(a) 

AB = 



2 

0 




-2 




2 

0 



, and C 





2 k 




-2 2 






0 1 




0 0 










— J 


t 





3 

1 



-4 4 

0 0 



(b) 



” -4 4 ~ 


■ 


* 1/2 3 




r ~ 


0 0 




-3 1 










- 




— . - 



(AB)C = 






I 



i 






& 



I 






V 



g 



£ 

1 






O 
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AB(C) = A(BC) is also true in general. It is similar to the associative 
law of multiplication of ordinary algebra where (ab)c = a (be). 



A * 


‘ 2 4 

0 1 


1 I 

, B = 


r- 

-2 2 

0 0 


, and C = 


1/2 3 

-3 1 


then 

(a) 


* 




k»— 







B *!- C = 



■2 

0 



2 

0 



1/2 

-3 



3 

1 



-3/2 

-3 



5 

1 



(b) 



A(B J r C) = 



2 4 




-3/2 5 




r 


0 1 




-3 1 






u. — J 


1 


— • 




■ mm 



-15 

-3 



14 

1 



(c) 


2 


4 


f -2 


2 : 


f 


— j 


AB = 


0 


1 


u 

1 


0 


i 

' L»... 











COMMENT 









i 



$ 

$ 



If A, B and C are matrices of the same size then A(B -I- C) = AB AC which 
is called left distributive law of multiplication with respect to addi- 
tion and is similar to the ordinary algebra distributive law (a(b + c) = 
ab -I- ac). 

And 



(B -I- C)A = BA v CA which is called right distributive law of multiplication. 
Since multiplication of matrices is not commutative, we cannot conclude 
that the left distributive multiplication is equal to the right distributive 
multiplication, i # e. , A(B -!- C) •■!■ (B -!• C)A, though it is true in ordinary 
algebra where multiplication is commutative and therefore 
a(b c) = (b c)a. 






1 



l 



I 

% 

I 
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11. Let 



A 



then 
(a) A(B 



m 



m 



mm 
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! 



2 

1 



-1 

3 




and C " 



2 





• • • • c 





and 

(b) AB 







o 

ERIC 









lalii 




This illustrates both the failure of the commutative law of multiplication 
and the breakdown of the law of cancellation. 
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I 



13. Let 





a b 




c d J 


A = 




and B = 






-b a 




-d c 



then 

(a) 



AB = 



a 

-b 



-d 



ac - bd 


ad + 


be 


be - ad 


-bd + 


ac 



and 

(b) 




ca - db 


cb + da 


-da - cb 


-db -!■ ca 



(c) We figured out that 









.. ... 








' v ‘ 




ac - bd 


a d -!* 


be 




. ca - db 


cb + 


da 


AB = 


-be - ad 


-bd -!• 


ac 


, and BA = 


-da - cb 


-db -h 


ca 



















It shows that AB and BA are 

(aqual/tmequal) 









100 



and 

(b) 





0 

a 




e • • • • 







(c) Therefore AB - -BA 

If AB = -BA, A and B are said to be anticommutative. 



16. Let 



A 




1 

3 



then 

AA - 5A + 21 

















2 1 




2 1 




2 1 




1 0 








-5 




-1-2 




4 3 




4 3 




4 3 




0 1 






- 






■ 


- ~ 





(name of the matrix) 












- 




r 




8 5 


+ 


-10 


-5 


» 

I 


2 


0 


20 13 




-20 


-15 




0 


2 


- -j 








t 




— 



0 0 



0 Oj 

It shows that A satisfied the equation AA-5A -!- 21= 0 
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17. Let 



~2 0 7 ~ 




-x -14x 7x 




'" 1 o o' 


0 10 




0 10 


=s 


0 10 


12 1 




x 4x -2x 




0 0 1 



then 







- 

10 0 






0 10 






0 0 1 


(multiplication of two 


ma 


- . \ • 
trices) 



5x 

0 



0 

1 



0 

0 



0 -10x‘i-2 5x 



13. If 



— * 

5x 


0 


0 




1 


0 


o’ 


0 


1 


0 


= 


0 


1 


0 


0 


- 10 x -:-2 


5x 




0 


0 


1 



then x = 



1/5 



\M* 7L 



sa&Baaai&aasaaa^^ 
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Unit IX: Powers of Matrices 



Since matrix multiplication and ordinary multiplication are similar in 



many important ways, we can define the powers of a matrix in the ordinary 
way . Our definition shall lead us to the basic laws of exponents . 

We define A 2 = A-A, A 3 = A’A 2 , A 4 = A -A 3 , etc. That is, we make the induc- 
tive definition: A 2 = A*A, = A*A n . 



It is also convenient to take A 1 = A and A° = I. 



1. Let A be any matrix. Then 



A 2 * A 3 = (A* A) (A ‘A-A) = A’A'A’A’A = A x 



What does x equal? 



x = 5 



2. Let B be any matrix. Then 

( B 2)3 « ( B .B) 3 = (B-B) (B* B) (B-B) 



= B‘B“B # B*B*B = By 



What does y equal? 



y = 6 



CJ 
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Let A be a matrix. Then 



A n -A m = A“’ 
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COMMENT 

It should be noted that the law of exponents applies only to powers of a 
single matrix A. We cannot conclude that 

0 0-3 

ABAB = A B , 

since this would involve rearrangement of factors. Such a formula can 
very well be false. Care must be taken to keep the factors in a product 
in correct order. 



I 



f 

I 

I 

i? 

1 



I 
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i 

3 
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I 

I 
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! 

>1 

I 
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On the other hand, a formula like 
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o 

O . 



Let 




0 

1 



0 

0 



0 0 0 



then 


10 0 


1 


’l 0 o' 


A 2 = 


0 10 




0 10 




0 0 0 




0 0 0 



(name of 
the matrix) 



10 0 
0 10 
0 0 0 



A 



2 

If a 



- B, 



B is said to be an idempotent matrix. 



Therefore in the above 



frame A is an idempotent matrix because A 



2 



= A 






9. If 




‘ 2 


-2 


-4~ 




2 


-2 -4* 


A - 




-1 


3 


4 


2 

then A = 


-1 


3 4 






1 


-2 


-3 




1 

L, 


-2 -3 


Therefore 


A 


is 


called 








.matrix 



idempotent 



COM M ENT 



Let A and B be matrices of the same. size. Then 
and 



(2A) (3A) = 2((3A)A) = 2(3(A*A)) = 6A‘ 
(2A) (3B) = 2((3A).B)= 2(3(AB)) - 6AB 
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Note that ve continue to pay careful attention to the order in which the 
matrices come, 6AB must be distinguished from 62A; only the numerical 
factors, and not the matrix factors, of a product may be rearranged. 








(6A)(3A) (B)(4) = 

(after multiplication) 
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Let A and B be two matrices of the same size and these do not commute, i.e. , 
AB 4 BA. Then 

(a: b) 3 » (a-: b) (a-i-b) (a-:-b) 

= (A*: B) (A (A-:-B)-:-B (A+B ) ) 

- (A-:B)( ) 

= A (A 2 -!-ABrI-BA-:B 2 )+B (A 2 -: AB-'.-BAH-B 2 ) 

= A 3 a 2 b aba ab 2 •: ba 2 v bab *:• b 2 a b 3 



i 



i 



2 2 
a ab-: ba+b z 



I 




f COMMENT 



A 



% 



We have seen that in matrix algebra 

(A: b) 3 « a 3 -ia 2 b-i-aba+ab 2 -:-ba 2 *: bab*:-b 2 a>; b 3 

And there we must stop. If A and B were numbers we could add A 2 B*i-ABA-:~BA 2 
to get 3A 2 B, and AB 2 -!-BABvB 2 A to get 3AB 2 . But this involves a rearrange- 
ment of factors which is permissible for numbers, impermissible for matrices. 
This simple counter-example rules out the possibility of there being a 
binomial theorem for (non- commutative) matrices. Therefore the break- 
down of the commutative law for multiplication for matrix algebra makes 
another difference between matrix algebra and ordinary algebra. In 
ordinary algebra where binomial theorem holds good (a*!b) would have been 
expanded very easily without multiplication, with the help of the binomial 
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theorem as follows: 
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(a-Hb) n = a n -:-na n " 1 b'l iiliidLl a n " 2 b 2 -!- -:-b n 

21 

(where ni = n(n-l)(n-2) 3 . . . . 3*2*1 

Using this binomial theorem 
(a-!-b) 3 = a 3 + 3a 2 b-K3ab 2 -i'-b 3 



\ 



14. Let A and B be two non-commutative matrices of the same size. Then 

(A+3B) 2 (A-2B) = ( ) (A-2B) 

. = (A 2j - 3AB-;-3BA-!-9B 2 )A-!- (A 2 *i-3ABr!-3 BA-1-9 B 2 ) (- 2B) 



a 2 -i-3Ab*:*3ba+9b 2 

A 3 -:-3ABA-:-3BA 2 -: 9B 2 A-2A 2 B-6AB 2 -oBAB-1GB 3 



15. Let A and B be two non-commutative matrices of the same size. Then 
(A+B)(2A+3B)(A - 3B) = (A-:-B)(2 A(a-3B)-:3B(A-3B)) 

= (A+B) ( - -5* - ) 

= A(2A 2 -6AB-!-3BA-9B 2 )-IB (2A 2 -6AB-1-3BA- SB 2 ) 



2A 2 -6ABv3BA-9B 2 

2A 3 -6A 2 B-!-3ABA- 9AB 2 H-2BA 2 -6BAB*}-3B 2 A-9B 3 
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16. Let 



1 



A = 



-1 

L. 



then 



-1 

1 



r— 



A 2 -2A - 



1 

-1 



.l] 



1 

-1 



• 1 
1 



-2 



1 

-1 



-1 

1 



r* 



(name of the matrix) 



a 



I 

1 



% 

% 



3 



1 

1 

I 



3 : 



$ 

4 ? 

% 



-2 

0 

0 



•2 

2 

0 

0 



-2 

2 



2 

-2 



0 



It shows that A satisfies the equation A -2A = 0 



17. Let 



A 



1 

-1 



0 

2 



then 



A 2 -3A-:-2I = 



1 

•1 



1 

■1 



0 

2 



-3 



1 

-1 



0 

2 



• 1-2 



1 

0 



0 

1 




:::::] 



(name of the matrix) 






SH 









i 

| 

I 

I 



$ 

a 



l 

Si 



I 



■>i 



i 



1 



| 

I 



j 








18. We saw in the previous frame that 

A 2 -3A -!- 21 « 0 ■* 

therefore A the equation A 2 -3A 21 = 0. 

(satisfies/does not satisfy) 



satisfies 



19. Let 



A 



then 



1 

0 

1 



0 

1 

1 



1 

0 

0 





1 


0 


r 


k 

1 


1 


0 


1* 


A 2 = 


0 


1 


oj 


5 

1 , 
1 


0 


1 


0 




1 


1 


0 




1 


1 


0 



(a) 



1 



% 

i 




l 



$ 
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5 

4/x 

0 

0 



0 





|— imf 






4x 




X 

<• 

i 

CM 

1 




-3 0 




J 

1 




J r 




5 




CM 

1 

X 

t 




0 -3 




L. 




L. _ 



= 0 



0 



2 

(b) Therefore, A -2A-3I = 0. 

It shows that A satisfies the equation 



A 2 -2A-3I = 0 

No matter what the value of x the matrix A satisfies the equation 
A^-2A-3I = 0, and therefore this quadratic equation has infinitely 
many distinct 2x2 matrices as roots. 

. * 



f 



h 



COMMENT 



We have seen that the law of exponents is the same for powers of a matrix 



as for powers of a number and that numerical factors in a product may be 
rearranged at will. Moreover, we know that the laws governing addition 



and subtraction of matrices are just like the laws governing addition 
and subtraction of numbers. It follows from all this that the algebraic 



rules for manipulating expressions which are made up out of sums of 



powers of a single matrix A, multiplied by arbitrary numerical coefficients , 






are just like the algebraic rules for manipulating ordinary polynomials. 



Thus, for instance, 
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(A-I-2I ) = A* A-i-A . 2I-1-2I . A-1-2I • 21 

= A 2 -!-4A-!-4I 



just as 



(x+2) 



= *2+ 



x -5-4x-!-4 



and 



(A-I) (A-i-I) = (A-I) (A -I) 



3 2 

= A -A -A-M 



just as 



(x-1) (x-I-l) * (x-l)(x -1) 



3 2 .i 

= X -x -x-.'-l 



21. Factor the following polynomials in A into a product of first-degree 



polynomials in A: 



(a) A -I = 



(A-I) (A-I-I) 



(b) A 3 -! 3A 2 -; 3A J rI o 



(A-I-I)' 



(c) A -4A-I-4I = 












Unit X: Determinants 



$ 



So far you have learned addition, subtraction, and multiplication 
of matrices. Now before proceeding to division, we will discuss an 
operation applicable to the elements of a square matrix that leads to 
a scalar value known as the DETERMINANT of the matrix. 

A determinant is a polynomial of the elements of a square matrix. 
The notation for the determinant of the matrix A is |A|. Determinants 
are defined only for square matrices — the determinant of a non-square 
matrix is undefined and does not exist. The determinant of a matrix is 
a scalar. Obtaining the value of j A} by adding the appropriate products 
of the elements of A is referred to as evaluating the determinant, ex- 
panding the determinant or reducing the determinant. 

1. The determinant of a 1 x 1 matrix 

A-0] 

is written as 

= M ' 



|Aj 

Note that the determinant of a lxl matrix is the value of its sole element. 



| 

I 

h 



I 



1 



I 







•-k 

1 



i 

■r. 

1 




2. The determinant of a 2 x 2 matrix 



B = 



11 b 12 



I 



21 



22 



_ 



is written as 



|B| = 






• • • *+ • 






erJc 



i 
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11 



12 



b 21 b 22 



3. The determinant of a 2 x 2 matrix consists of the product of the diagonal 
terms minus the product of the off-diagonal terms. Hence in general 



M = 



b 

d 



3(3 “ • • • • • 






be 



4. Tlie determinant of a 2 x 2 matrix 



M = 



3 7 

17 20 



••••• “ •••*. — 59 



3x20 - 7x17 



5. The determinant of a 2 x 2 matrix 



10 1.6 




-92 -8 = -100 
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6. Evaluating the determinant of a 2 x 2 matrix- is known as the evaluation of 









8. In the previous frame second-order determinants 



a 22 


a 23 


• 

9 


a 21 


a 23 


; and 


a 21 


a 22 


a 32 


a 33 




a 31 


a 33 




a 31 


a 32 



which are obtained from the determinant j A | , are known as minors of the 
elements , a ^ and a^g respectively, where a.^, an< ^ a 13 are 

the elements of the 1st row of the determinant |Aj. 



9. The evaluation of the determinant of a 3 x 3 matrix 






3 1 




2 1 






II 

/-S 




+5(-l) 




-:-6 (-: 1) 






9 7* 1 




CO 







= 4( - ) -5(. . )-:-6( - ) 



= 4xl2-5x6-I-6x-.. . . .=-18 



2 3 

S 9 

(21-9); (14-8); (18-24) 

-6 

10. In the previous frame: 

(a) Second-order determinant 
3 1 

9 7 



is the minor of the element 



in the determinant j A 







I 



5 6 




4 6 




4 5 




• 

> 




• 

9 




3 1 




2 1 




2 3 



( 5 - 18 ); ( 4 - 12 ); ( 12 - 10 ) 



COMMENT 

The method of expanding a determinant which we have learned in the 
previous frames by using elements of any row can also be applied to the 
elements of any column. No matter by what row or column the expansion 
is made^ the value of the determinant is the same. This method of ex- 
panding a determinant is known as expansion by the elements of a row 
(or column ) or as expansion by minors . The (-M) and (-1) factors are 
decided on according to the following rule: if A is written in the form 
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i : 



124 




1 



£ 



s 



I 



K 



£ 



P 



ERIC 









A = fa 1 , the product of a., and its minor in the expansion of the 

t i jJ 1J 

^ i-:-j 



determinant | A | is multiplied by (-1) • Once a row or column is 

decided on and the sign calculated for the product of the first element 
therein with its minor, the signs for the following products alternate 
from plus to minus and minus to plus. The product of f.**l) to the 



corresponding minor is known as the cofactor of in |a|. 



13. The expansion of a fourth-order determinant is an extension of the ex- 
pansion of a third-order determinant which you have learned. Thus in 
the expansion of a fourth-order determinant 



A I - 



1 


2 


3 


4 


5 


6 


7 


8 


9 


10 ' 


11 


12 


13 


14 


15 


16 



by the first row 

(a) the minor of the element 4 is 



5 


6 


7 


9 


10 


11 


13 


14 


15 
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(b) The product of the element 4 with its minor is multiplied by 

(-I)* • • # __i 



(-D 



1:4 



14. The expansion of the determinant 



l A l 



1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 



by the second column is 



} A j = 2 (- 1 ) 


• • • • • • • • • 

9 • • ••• ••• 




• • • ••• ••• 

• •• ••• ••• 


-'-10 (- 1 ) 


• •• ••• 

• •• ••• ••• 




• « # • • # ••• 




• m • • • • ••• 




. • • • ••• ••• 



l4(-;-l) 



= -2i5x! 



I 



« j , { 


• • • • • * 




• • • • • • . 


| 


• • • • • • 




• • • • • • 


-7x 




+8x 




>-!-6< lx 

4 \ 




-3x 




T • • • •••) 


• • • • • • 




• • • • • • 


t 


• • • • • • 




• • • • • • 



‘1-4 x 



i r . . 

ft ft • • • • • • 1 

r 10'lx -3x 


• • • * • • I 


1 

-1-4 x 


• • • • • • 


LL 


• • • • • • 


1 1 • • • • • • I 

J v 1 


• • • • • • 




• • • • • • 


[ l 


• • • • • • 



*v 

| • • • ••• 5 • • • t • « I I 

-3d -;4d I K 

I • • • ••• I*.. 9 • • I jl 



=-10(176-180)-’.T4(144~156)-16(135-143)-!-6 (176- 180) -13 (144- 156) 
24 (135- 143) -10 (112- 120)-!-30 (80- 104) -40 (75-91)-! 14(84-88) 
-42(60-72)-!-56(55-63)=0 









1 



3 



4 



i 



I 



*■ 



i 



ERIC 






2<-D* 



5 7 8 

9 11 12 

13 15 16 



•!-6 (•:•■! ) 



9 11 12 

13 15 16 





13 4 




1 3 4 


-:-10(-l) 


5 7 8 


-:-14(+l) 


5 7 8 




13 15 16 




9 11 12 



r ! 


j 11 12 




9 12 




9 11 


1 f 


11 12 




9 12 




9 11 




\5x 




-7x 




vox 




?+6l lx 




-3x 




-i-4x 




l 


l ! 


1 15 16 




13 16 




13 15 


l 


15 16 




13 16 




13 15 





-lOx! 



7 


8 


-3x 


5 


n 

O 


-l-4x 


5 7 


> +14 jlx 


7 


O 

O 


-3x 


15 


16 




13 


16 




13 15 


i 


11 


12 





5 

9 



O 

o 



12 



-!'4x 



5 7 
9 11 



COMMENT 

The method of expanding a determinant used in the previous frames 
for expanding determinants of third-order and fourth-order, can be used 
for the expansion of determinants of any size. The determinant of an 
n x n matrix A = [ a ijj °htained, hy considering the elements of any 



one row (or column), as follows: multiply each element, aij , of any 

row (or column) by its minor, |m^| > t^e determinant derived fromjA j 



by erasing out the row and column containing a^.. * mu ltipiy the product 



by (-1)^'^; add the signed products and their sum is the determinant| A 
This expansion is used recurrently when n is large, i.e. , each |m^| 
is expanded by the same procedure. 



15. Let 



A = 



■ 3 

-1 



6 



2 

1 

7 



4 

5 

10 



; then 



(a) 



✓ 

A = 






128 



\ 








—2 4 




7 4 




7 -2 


= 3x 




-lx 




-f-5x 






9 10 




6 10 




6 9 



= 3(-20-36)-(70-24)-:-5 (63+12) 



00 




7x 


• • t • 

• • • • 


• ♦ w • 

• • • • 


-1-2 x 


• • • • • 

J 

• ••• • • • • 


: 

-!-4x 


7 (. . 


- 


.)+2(.. 




>I4( 


) = - 161 



• • • • 



• • • • 





1 5 




3 5 




3 1 


= 7x 




-:-2x 




_;.4x 






9 10 




6 10 




6 9 



= 7(10-45 >:-2(30-30)-J-4(27-6) 



(c) From (a) 

I A | = 161 
ana from (b) 

j B | = -161 

therefore j A j = ( ) | B | 



( - ) 
minus 

Note : Notice that by changing two rows of the matrix A, the sign of its 
determinant is changed. It is true in general that interchanging 
two rows of a determinant changes its sign. 
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17 . The 


expansion 


of 








a ll 


a 12 


a 13 


M 


> 


a U 


a 12 


a 13 






g 31 


a 32 


a 33 



is | A | 




a 12 a 13 




a ll 


a 13 




a ll a 12 


= a ll 


a 32 a 33 


” 3 12 


a 31 


a 33 


-J*a , 
13 


a 31 a 32 





a ll( “ )” a 12( “ )" ;_a 13 ^ " ) 



JL 

k ■ * 



1 

k * 4 



a ll (a 12 a 33" a 13 a 32^“ a 12^ a ll a 33“ a 13 a 31^ : ' a 13 (a ll a 32“ a 12 a 31 ) 

a ll a 12 a 33“ a ll a 13 a 32" a 12 a ll a 33 +a 12 a 13 a 31 H ' a 13 a ll a 32" a 13 a 12 a 31 
= 0 



Note : Notice that here | A | = 0. It can also be seen in another way: 

if the first two rows of |a| are interchanged then the value 
of the determinant is unaltered but the sign is changed, that is, 
| A| = - j A j so that | A| =0. It is true for general that if two 
rov/s (or columns) of a determinant are the same the determinant 
is zero. 



18. Let 




12 18 



A 
/ r* 



1 




then 



6 



7 



9 



(a) 

|A| 





3 


12 


18 




1 


5 




2 


5 




2 1 


s 


2 


1 


5 


= 3x 


7 


9 


-12x 


6 


9 


-!-13x 


6 7 




6 


7 


9 


= 3 ( t . 




.)-12(... 






>: 18( — 





= 210 
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3(9-35)-12(18-30)-;-MI4-o) 
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I 
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£ 



j* 
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(b) also 





1 4 6 


r 

i 






l&J = (3) 


2 15 


=3 j lx 


i 

1 


-4x 




6 7 9 


L 


• • • • «* • ] 





= 3 £(9-35)-4(lG~30)-:-6(14-6)J =210 






1 5 




2 5 




2 1 


lx 




-4x 




-1-6 x 






7 9 




6 9 




6 7 



Note: If any scalar is a factor of a row (or column) it is also a factor 

of the determinant. As in the above case 3 was the factor of the 
first row and it is also the factor of the determinant. 



19. Let 


a ll a 12 a 13 




a ll a 12 a 13 


IM - 


3a ll 3a 12 3a 13 


I! 

W 

sy 


Q 11 a 12 3 13 




a 31 a 32 a 33 




a 31 a 32 a 33 



Therefore the value of I A | =0, because the first two rows of the right 
hand side determinant are 



Note: 



the same or equal 

It is true in general that if one row (or column) of a determinant 
is a multiple of another row, the determinant is zero, as factoring 
out the multiple reduces the determinant to having two rows (or 
columns) which are the same and hence it is zero. 
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23. Let A be an n x n matrix and y is a scalar which is multiplied to each element 
of A, then | yAj = y # * * * j Aj 



n 



24. Let 



i a| = 



2 

-4 

1 



3 

4 
7 



= 1 (-28-4) -2(35-36)+3 (5+36) * 93 



then multiplying the first row by 3 and adding it to the third row 



we get 





1 


2 


3 




1 


r\ 

£ 


3 


|A| - 




-4 


4 


S 


5 


-4 


4 




9+. ... 


1+. ... 


7 +. ... 











= l(-64-23) -2(00-48)+3(35+4d) = 93 



1 


2 


3 




1 


2 


3 


5 


-4 


4 


S 


5 


-4 


4 


9+3 


1+6 


7+9 




12 


7 


16 



Note: 



Notice that when three times the first row is added to the 3rd row 
the value of j Aj remains the same. It is true in general that adding 
to one row (or column) of a determinant any multiple of another row I 
(or column) does not affect the value of the determinant. 
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25. 



We have seen that adding a multiple of a row (or column) to another 
row (or column) does not affect the value of a determinant, but adding 
a row (or column) to a multiple of another row (or column) is not the 



same thing a 


nd 1 


eads 


to a 

j 


different result. 






a ll 


a 12 


3 13 






1 A l 


1 = 


a 21 


a 22 


a 23 










a 31 


a 32 


a 33 






then 


adding row 2 to 


y times row 1 we get 




y a n* : 


a 21 


• • • 


J- 

1 • • • 


-X. 






a 21 




a 


22 


3 23 






a 31 




a 


32 


a 33 




which 


; is equal 


to y 


i*i 


which is different 



y a 12*' a 22 » y a 13 +a 23 



C 0 M M E N T 

So far we have learned the following properties of determinants: 

1. The determinant of the transpose of a matrix is the same as the 

determinant o/f the matrix: | A j = | A | 

2. If two rows (pr columns) of a matrix are interchanged, the determinant 
of the matrix changes sign. 

3. If two rows (or columns) of a matrix A are identical, then I A | = 0. 

4. If all the elements of a row (or of a column) of a square, matrix A be 
multiplied by an element k, the determinant of the matrix is multiplied 
by lc. 
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26. 
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5. If two rows (or columns) of a matrix A are proportional, then | A | - 0. 

6. If a matrix A has a row (or column) of zeros the value of its deter- 



minant is zero: 



A 



= 0x| A |. 



7. 



6. 



If A is an n x n matrix and lc is a scalar then jlcAj = k n jA|. 

Adding to one row (or column) of a determinant any multiple of another 
row (or column) does not affect the value of the determinant. 

These properties can be applied in endless variation in expanding 
determinants and thus the expansion becomes much simpler. 



The expansion of the determinant 
15 3 



i*i ■ 



2 11 7 

4 21 15 



= 1 (165-147 )-5 (30-28)-: 3(42-44) » 2. 



| A | can also be expanded as follows: 



If -5 times column 1 
-3 times column 1 



is added to column 2 
is added to column 3 



and 

we get 



I a| - 



1 0 



2 1 



4 1 



« lx 



1 1 

1 3 



= 2 



0 

1 

3 
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Unit XI: Inversion of Matrices 



In the previous units of this program we have studied addition, subtraction, 
and multiplication of matrices and the properties of determinants. It 
shoul.d not be surprising that our next aim is to study the division of 
matrices. In ordinary algebra, division is necessary to solve equations 
like ax=b. There is a similar problem with matrices, viz., solving equa- 
tions like AX=B. Since the commutative law of multiplication breaks down 
for matrices, AX need XA. 





not equal 


1 


problems. 


different 




Y 



not equal 



mmmmm 



i Sifna. jycri aarawnaitesir r. z kz* im r** ** -. i .- rt . 41^*1**"' 
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r 








— * 


7. Now let 

A = 


1 C 

1 1 


, 3 = 


, 

1 2 
0 1 


, and X = 

* 


-1 2: 

'I 1 




A 




— — 




- J 



then 
XA = 



r 




— ^ 




-1 2 




1 0 


r 


-1 1 




1 1 


1 •••••_ 


- 




~ — 





(name of the matrix) 



r 

1 



o 



= B 



Notice that X = 



- 1 
— * 



-1 



2 

1 



is the solution, of the equation XA-B. 



8* Mow we have seen that the solution of the equation 

I 1 2 



AX— B ' is j 



V 

1-1 



XA=B is 



-1 



and the solution of the equation 



It shows that the value of X in both cases 



is 



( s ame / differ enc ) 



different (not same) 



Notice that it is due to the fact that AX /' XA. 
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Notice that the value of X = 



!-l 



11 

-9 



3 1 

-2 -1 



G 0 M H E N T 

Nox 7 we have seen that if we multiply through (on the left) the equation 



AX = B by another matrix A"^' such that A a = I, then we can solve the 



equation as follows: 



-1 



AX 


= B 


AX 


- A'*: 


IX 


- A' 1 


xr 

A 


= a" 1 : 



Thus X = A~ is the solution to the equation AX = B 

Since behaves very much like a reciprocal in ordinary algebra, 
we shall call it the reciprocal matrix of A or the inverse matrix of A, 
Thus the problem of division of matrices is actually the problem 
of finding the reciprocal or inverse of a given matrix. Some matrices 



have reciprocal matrices, some do not. 



To find out the inverse of a 2 x 2 we proceed as follows: 



r 

? 

;a 
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1° 



I 



and its inverse A 



-1 _ 



jP 

r 



q 

s 



Then 



i 


r m 

1 




f i 


a b ' 

i 


t 


p q 




i o! 


c d 




r s 




0 1 


- 




L 







if A has an inverse 



or 



| ap br 
i cp + dr 



aq + bs 
cq + ds 



i0 

L 



o! 

l 






mLmmm 




.* .. - ,. , .„, ' .... , 4 - - . 
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15. If the inverse of a matrix A = 



b 

d 



is 



equal A"^ = 



-b 



ad -be ad -be 



-c 



ad -be ad -be 



Then 



(a) The inverse of the matrix B = 



2 

1 



4 

4 



is 



S' 1 is 



• * • • • 



r~ 



-1/4 



-1 

1/2 



(b) The inverse of the matrix C 



0 

1 



***> 



2 



is 



_ i 

C x = 



-2/3 1 



1/3 



0 
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18. Let 



A = 



0 

1 



3 

2 



, A 



-1 _ 



Solve the equation AX = B 



-2/3 

1/3 



1 

0 



and B = 



1 

-1 



1 

2 



i 

i 

I 











-2/3 1 




i 

1 1 




-5/3 4/3 


1/3 0 




-1 2 




1/3 1/3 


— mm 











19. Let 
A 



1 

0 



0 

0 



0 

1 



0 



Then 



1 0 





1 


0 


0 




1 


0 


0 




r ~ 


AA = 


0 


0 


1 


1 


0 


0 


1 


= 






0 


1 


0 




0 


1 


0 







(name 
of the 
matrix) 



1 0 
0 1 
0 0 









Notice that AA = I, therefore A is its o\m inverse 



20. Let 



1 


0 


0 


|C •' 


2 . 


1 


0 


0 


0 


1 


~ and B = 


1 


2 


0 


0 


1 


0 




2 


1 


1 



Then 

Solve the equation AX = B 




1 


0 


0 




2 


1 


-M 

0 




2 


1 


0 


0 

L° 


0 


1 




1 


2 


0 


— 


2 


1 


1 


1 


0 




2 


1 


1 




1 


2 


0 

«4M| 



21. Let 



A 



Then 



1 


0 


1~ 




1 


0 


-1 


1 


1 


0 


and B = 


-1 


1 

X 


1 


L° 


0 


1 


i 


0 


0 


1 



(a) 


1 0 1 ~ 




10-1 


1 


m- 


AB = 


110 




-111 


= 






0 0 1. 




I — 
o 

O 

l 







(name 
of the 
matrix) 
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1 

0 

0 



0 

1 

0 



0 

0 

1 



Notice that when AB = 



I, then A is the inverse of B, or A = B”^ 






(b) and 


10-1 




1 

t— 4 
O 

l — 1 

... i 






BA = 


-11 1 


\ 


110 


s 






0 0 1 




0 0 1 








* — ■ 1 


Lm ^ 


i 





(name 
of the 
matrix 



i 

0 

0 



0 

1 

0 



0 

0 

1 



Here BA = I, therefore B is the inverse of A, or B = A 



-1 



C 0 M M ENT 

T-Te have seen in the previous frame that AB = I = BA which shows that A 
is the inverse of B and B is the inverse of A. if we put B - A"-*- in the 
above equation we get 



A A-l = I = A _1 A 



Therefore A commutes with its inverse. It is true in general that every 
matrix commutes with its inverse provided the inverse exists and is similar 
to the ordinary algebra where aa***- = 1 = a“^a. Also if a matrix has an 
inverse, it has only one inverse; that is, this inverse is unique. 















3 




3 



5 




1* 



l 






V 

£ 
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Then 

(a) A 2 




2 

5 




2 

5 







% 



| 



I- 

I 



r* 



I 13 

t' *** 

i 



(b) 


— — . 




■» — j 








1 2 




1 0 




• •••• ••••• 


A - I = 




- 




rs 






3 5 




0 1 








_ .. 




L J 


^ — 



12 

31 



I 



I 

I 

I 

I 

i' 




l 



0 2 

3 4 



(c) 

A 2 -7I 









•r** 










7 


12 




1 


r\ 

u 








i 


-7 




= 






10 


31 




0 

'4M* 


1 








%/ X 








- 





0 



12 



I O 

o 
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(d) 








wi - 




i 


0 


12 




0 


2 


(A -7t))-S(A-I) = 


18 


24 


-6 


3 


4 


1 








— 


— 



:]■ 



23. We saw in the previous frame that 
(A 2 -7 1) -6 (A- 1) = 0 



0 

0 







$!•' 

? 

K 






or A^-6A-I = 0 

Therefore the matrix A satisfies the equation 

= 0 



% 



I 



I 



$ 



E 



A2-6A-I = 0 



C 0 M M ENT 



Once we have discovered an equation which is satisfied by a particular 
matrix, it can be used to find a reciprocal for the matrix. As the equation 
A2-6A-I = 0 can be written as 



A2-6A — X 



or A(A-6I) = I 



Then (A -6 I) must be the inverse of A, 
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24. If A = 



1 

3 



A 



-1 



2 
5 

- I 1 



and (A-6I) is its inverse, then 



A-6I = j 3 



-6 



0 



1 

-L 



& 







i? 



I 

I 



I 

| 





i 
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25. If 



A = 



Then 








I 



I 



1 

0 



0 

1 



26, 



Let A = 









1 0 


n 


1 ° 




, then A*" = 




1 l 




2 I 


- - 




L J 



and 







■9— 






— 








1 


0 




1 


0 






A-I = 


- 




- 






= 






1 


1 




0 


1 










1 


1 




— 
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-A-1-2I 



28. 



If A 2 -2A-rI = 0, then 
A 2 = 2A-I 
4 



So A' = (2A-I) • • • 

= 4A 2 - 4A + I 



4A2-4A+I can further be reduced by substituting 
A 2 - 2A-I, getting 



4(2A-I) -4A-1-I 



« 4A-3I 



A 



2S. If A - 4A-3I, then 



A S = (4A-3I) 2 



16A 2 -24A-!*9I 



30. Polynomial . 16A 2 -24A+9 I can be simplified further by substituting A 2 = 2A-I, 
getting ; 



16(2A-I) -24A-I-9I 
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8A-7I 



C C H M ENT 



The polynomial equation satisfied by A is useful in doing other calcula- 



tions involving A. As we obtained the value of A° without much calculation. 



Once the polynomial equation of a matrix is obtained, it is much easier to 
study all **he other properties of the matrix. 

Let A be an n x n matrix. Then A satisfies a polynomial equation of 
the form 



A k -!-c k _ 1 A k ’ 1 -i- -!~c 0 I = 0 



of a degree lc which is not more than n and with numerical coefficients 



'fc-1 » 



> C Q. 



Thus, a 2 x 2 matrix always satisfies either a linear or a quadratic 
polynomial equation . A 3 x 3 matrix always satisfies a cubic, quadratic, 
or linear equation . A 4 x 4 matrix always satisfies a polynomial equation 



of degree 4 at most, etc . 



Let 




■4*1 
















' A = 


1 


3 


, then A“ 






4 


6 








1 


1 








2 


4 






and 
















- 




— — ^ 


(a) A- I = 






- 




= 








- 










- 



3 

1 



1 

0 



0 

1 



0 

1 



3 

0 
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(b) Therefore 

A 2 ~9I 









9 


-6 


-9 


i 

X 


A 

U 


-30 


69 




0 


1 




* 


M. mJ 



-30 60 



Wow 


•*» 








-*» 


A 


0 


-6 




0 


1 


A -91 - 


-30 


60 


and A- 2 1 = 


c 


-10 






•« 






* 



(a) Fhat multiple of (A-2I) should be subtracted from or added to 
(A 2 -9I) so that all the elements of the resultant matrix may be zero? 



6 times added 



(b) Therefore 

A 2 -9 1+6 (A- 2 1) = 



0 

■30 



-6 

60 



-:-6 



0 1 
5 -10 



• ••••» 









, 


0 0 








= 0 




.0 Q 
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39. When A2 -91+6 (A -2 I) = 0 
or A2+6A-21I = 0 

Then the matrix A satisfies the equation 







A 2 +6A-21I = 0 



C 0 M WENT 

From the previous frames you must have got an idea of how to find an 
equation satisfied by a 2 x 2 matrix. To summarize it, let A be a 2 x 2 
matrix and I a unit matrix of size 2. To find the equation which may be 
satisfied by A, first we calculate A2. Then find a number "a such that 
A-al is a matrix with element 0 in the upper left-hand corner. If it 
happens that A-al is the zero matrix, we are done, for in this case A 
satisfied the linear polynomial equation A-al = 0. If, however, any 
element of A-al is different from 0, we continue the process by finding 
a number ! ’b : ' such that A 2 -bI is a matrix having 0 in the upper left-hand 
corner. We now have two matrices which have 0 in the upper left-hand 
corner; viz., A-al and A 2 -bI. Now, for 2x2 matrices, if A-al j- 0, 
then there exists a number !, c" such that 
(A2-bI) -c(A-al) = 0 

which is the required equation. 
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40. Let 

A = 



(a) To find a quadratic equation satisfied by the matrix A, first I 



should be multiplied by •••••••• and then subtracted from •••••••• giving 

(number) - (matrix) 



A-I = 



1, A, 



: 



(b) After I is multiplied by 1 and subtracted from A, the next step 
will be to multiply I by __ 

.... and then subtract it from where A 2 - 



(number) 



2 1 
1 1 



, giving 



A 2 - 21 = 



2 , 



0 

1 



-1 



(c) If A-I = 


'0 

1 


r 

-l 


and A2-2I - 


‘0 

1 


V 

-1 










m 





Then to find a quadratic equation satisfied by the matrix A, (A-I) should be 



multiplied by and then subtracted from (A 2 -2I) , giving 

(number) 



A 2 -2I-(A~I) = 



:] 
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1 , 



0 0 
0 0 



0 



(d) If A 2 -2I-(A-I) = 0 
or A 2 -A- I = 0 

Then this equation is satisfied by the matrix 



A 



41. If A 2 -A- I = 0 
or A 2 -A = I 

then the inverse of A is 



A-I 



42. If A 2 -A-I = 0 
or A 2 =A-KE 

then A^ : '= 

and after substituting A^= A-KE, it becomes 



A 2 -r2A+I 

3A+2I 
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43. If A 4 = 3A-1-2I 

then A 8 = 

and after substituting A / "— A+I, it becomes 
A 8 = 



I 



9A 2 -a2A-'r4I 

21A+13I 



r 



% 



l 



l 



v 



44. Now if 



A = 



1 

1 



and A° = 21A+13I 



mm 




Ji 






— 




l 


l 




1 


0 


= 21 


i 


0 


+13 


0 


1 
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13 



I 



45. Let 



A = 



1/3 

1/7 



(a) 



To find a quadratic equation satisfied by the matrix A, first I should 



be 



multiplied by and then subtracted from 



giving 



A- (1/3) I = 



] 
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(b) When I is multiplied by 1/3 and subtracted from A, the next step 
to find an equation satisfied by the matrix A will be to multiply I 



by 



o 

and then subtract it from A where 



A = 



giving A 2 - (10/9) I - 
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Then to find an equation satisfied by the matrix A, (A-(1/3)I) should e 
1 tip lied by and then subtracted from (A 2 -(10/9) I) , giving 
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A 2 - (10/9) I- (3.0/3) (A- (1/3) I) = 
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(d) If A2- (10/9) I- (10/3) (A- (1/3) I) 
or A'- (10/3) A = 0 
then A satisfied the equation 
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A -(10/3)A = 0 

Notice that A has no inverse as the equation which it satisfies has no 
term which contains I (alone) . 

1/3 7 



A = 



1/7 



has no inverse is true also if we test it from our earlier 



formula as follows: 



d/3) (3) -(1/7) (7) = 0 
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Unit XII: Inversion of Matrices (continued) 



In the previous unit you learned a procedure to find an equation 
satisfied by a 2 x 2 matrix A. This equation was used to find an inverse 
of the matrix A. The polynomial equation satisfied by A was also used to 
calculate A raised to a high power. 

The above procedure can be applied to matrices of higher orders. This 
procedure is based upon a theorem which tells us that every n x n matrix A 
satisfies a polynomial equation 



with numerical coefficients, the degree k of this equation being at most n. 

In this unit you will learn the procedure to find the least equation 
satisfied by a matrix A of any size and then find the inverse of A with the 
help of the least equation satisfied by A. 
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COMMENT 

We noticed in the previous frames that the matrix A of order 4 satisfied 
equation of 4th degree as well as equations of 3rd degree and 2nd degree. 
Therefore the equation of the second degree is the equation of the lowest 
degree which is satisfied by the matrix A. It should be clear that the 
equation of lowest degree which a matrix satisfies is the equation which 
gives the most information about the matrix. There is one and only one 
such equation of least degree for each matrix. The polynomial of smallest 
possible degree which a matrix A satisfies is called the least equation 
satisfied by A, 
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(b) 
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8, 
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If A * 
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(a) What multiple of I should be subtracted from A so that the element 
of the first row and first column of the resultant matrix may be zero? •< 
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(b) Therefore 
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(a) What multiple of I should be subtracted from A 3 so that the element 
of the first row and first column of the resultant matrix may be zero? • 
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11. Also 
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(a) What multiple of (A-I) should be subtracted from (A 3 -4I) so that the 
resulting matrix has zeros in row one column two as well as in row one 
column one? ••••••• 
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(a) What multiple of |(A 2 -2I) -(A-I)^ should be subtracted from 
f(A 3 -4I)-2(A-I^ so that the resulting matrix is the zero matrix? 
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(b) Therefore 
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13. We have 

|(A 3 -4I)-2(A-I)) - |(A 2 -2I)-(A-I)| * 0 

— ^ L ^ 

or A 3 -A^-A-I = 0, which is the least equation satisfied by 
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COMMENT 

In the previous unit we gave a systematic procedure for finding the 
least equation satisfied by a 2 x 2 matrix. 

Also in the previous frames you must have noticed that the same proce- 
dure was used for finding the least equations of 3 x 3 matrices. The 
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| preceding procedure, being general, can be applied to matrices of any size, 

2 3 n 

Let A be an n x n matrix. We first calculate A, A, A , We 

then find numbers a\ 9 a 2 , , an such that the matrices Bi - A-aj.1, 

B 2 = A^-a 2 I, , B n - A n -a n I all have the number 0 in the upper left- 

hand corner. If B^ =* 0, we are done, for then A-a^I = 0 is the least equa- 
tion we are seeking. If not, we find a second set of numbers b2, b^, 

b n such that the matrices C 2 - B 2 -b2®l> ^3 85 B 3“^3®1^ c n “ B n“^n B l 

have zeros for elements [ C] 1,2 as well ’as for [ C] x , 1* 

2 

If C 2 = 0, i.e, , if (A -a 2 I) -b 2 (A-a 1 I) - 0, we have the least equation satis- 
fied by A; if not the process is continued. 
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(a) To find the least equation satisfied by the matrix A, 
multiplied by and then subtracted from 
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(d) If A 2 + 10I-/2(A-/2I) = 0 
or A 2 -J2A + 121 = 0 

Then the least equation satisfied by A is 
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15. If A 2 -/2A + 121 * 0 
or (1/12) (/2I-A)A - I 
then the inverse of A is 
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A 2 -/ 2 A + 121 * 0 
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16. Let 
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17. Let 



A = 
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(a) To find the least equation satisfied by the matrix A , first I should be 
multiplied by •••••••••• and then subtracted from ••••••••••• } giving 
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(b) When I is multiplied by 1 and subtracted from A, the resulting equation 



A-I/0, therefore the next step to find the least equation satisfied by A 
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xtfill be to multiply I by •••••••••••••••••. and subtract it from A } 
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(c) Also multiply I by . 
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5 2 4 
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and subtract it from A 



where A J = 
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The next step to find the least equation is to get a second matrix (of the 
above two) with zero as element in row one, column two. Therefore (A -31) 
should be multiplied by and subtracted from (A 3 -10l) , giving 



(continued) 
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20. If 
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Therefore the inverse of A is 



A” 1 a 






Notice that A^ a j implies that A is its own inverse. 
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COMMENT 

We saw in the previous frames how to calculate the least equation satis- 
fied by a matrix. Once we know the least equation satisfied by a matrix A^ 
we can easily see whether or not A has a reciprocal. 
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If the final coefficient Cq in the least equation satisfied by A is different 
from zero, then A has a reciprocal as we can see from the following: 

We may write the least equation satisfied by A in the form 

+ CjA « -cqI 

/A .v * . .k-2 , . . 

or (A + c k _jA + + c^I)A = -CqI 

Since the numerical coefficient cq is different from zero, we can divide both 
sides of the equation by the number -Cq and find that 
-(l/c 0 ) (A k-1 + c k-1 A k " 2 + + cjI)A = I 



A k + c k . 1 A k ’ 1 + 
,lc-l 



This makes it plain that the matrix 

[-(l/ c 0> (A k 1 + c k _iA k 2 + + c 1 I) ] 

i 

is a reciprocal of A. 

If the coefficient cq in the least equation satisfied by A is equal to 
zero, then A does not have a reciprocal. 
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and let A“^ = B. To find the least equation satisfied by B, I should be mul- 
tiplied by •••••••• and added to B, giving 
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22, B + 51 jjfr 0, therefore to find the least equation satisfied by B, now I should 
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be multiplied by and subtracted from B" 
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23, We have 
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To find the least equation satisfied by B, (B + 51) should be multiplied by 



and added to B^-31I, giving 
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That is, the inverse of A is A, 
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We saw that if A and B both have inverses, then so does their product 
AB. Also it is true for general, i. e. , given a finite number of (square) 
matrices of the same size, each of which has an inverse, the inverse of the 
product of the matrices is equal to the product of their inverses written 



in the reverse order. 
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35. If A = 
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then AB 58 
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36. , | 
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ss b^A - 1 shows that 
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the inverse of (AB) is equal to the product of the inverse of 
and the inverse of 



B , A 



37. If (AB)~^ is the inverse of (AB) > then 
(AB)“ 1 (AB) » “ (AB) (AB) 
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If A satisfies the equation A 2 -2A+3I <= 0, does the inverse of A 

exist? 

(yes or no) 
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Unit XIII: Inversion of Matrices (other method) 

The inverse of a square matrix can also be derived from properties of 
the cofactors of elements of its determinant. A property of cofactors 
important for developing the inverse of a matrix is that the sum of products 
of the elements of one row (or column') with the cofactors of the elements 
of another row (or column) is zero. 



1. The cofactors of the elements of the first row of the determinant 
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(f) 


The sums 


of products 


of the elements 


of the 1st, second and 3rd row 


with 


the cofactors of the 


elements 


of the 


third row are: 
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1x76 + 4x-25 + (-8)x-3 = 0 
2x76 + 5x-25 + 9 x-3 ■ 0 
3x76 + 6x-25 + 10 x-3 * 48 



4. In frames 1, 2 and 3 we have seen that the sum of products of the elements of 
one row of A with the cofactors of the elements of another row of A is zero. 
Also the cofactors of the elements of A can be represented in a matrix B 
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(a) Each element of matrix B is the cofactor of the corresponding element 
of matrix A, Therefore 34 is the cofactor of 
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(b) The transpose of B is called the adjugate, or the adjoint, matrix of A« 
Therefore adjugate matrix of A is 
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(c) The inverse of A is calculated as 1 x B, 
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(d) 
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0 

0 

43 



« I 



Let 



A 



3 

6 

10 



1 

4 

5 



2 

5 

9 



(a) The cofactors of each element of A are: 



1st row 
2nd row 
3rd row 



11 , -4 



-10 



-3 



-3 



6 



b h 



-5 










% 







I 



l- 



l 



(c) The adjugate matrix of A » j 

L 



(d) and therefore the inverse of A & A"* 
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1 _ 

11 



-5 



3 



COMMENT 

Me have learned that if A is a square, nonsingular matrix its inverse, 
has the following properties. 

The inverse commutes with A, both products being the identity matrix: 
A" 1 A oI« AA" 1 . 

The inverse matrix is nonsingular. 

The inverse of A~* is A: (A~*)~* = A. 

The inverse of a transpose is the transpose of the inverse: 



(A)* 1 



(A* 1 ) . 



(5) The inverse of a product is the product of the inverses taken in reverse 
order, provided the inverse exist: (AB)“* ■ B~^A’*. 









i 



n 



4. 



A matrix of 2 rows and one column will be called a column of 

size 2* 






I 









% 



l 



vector 
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5* A matrix of size 3x1 will be called a column •••••*•••• of size 



vector, 3 



6. A matrix of size n x 1 will be called 



vector of size 



column, n 



7, The following matrix 



3 

2 

1 



is a column 



of size 



8* The following matrix 




is a 



of size 



vector, 3 



row vector, 3 







I. v - f 1 , 2, 3 "] is 



an example of a row vector of size 



10. v a LaJ is an example of a column vector of size 



11. Two row vectors , or two column vectors f are said to be equal if and only if 

corresponding components of the vector are equal. Thus the vectors 

v » [3, 4 I and u a r3, 4*1 are 

~ *■* (equal/unequal) 



12 . 





"l " 




»* •* 

1 


If U a 




and v a 






2 




2 




»3 . 




- 3 - 



then the vectors u and v are 



equal 



(equal/unequal) 



equal 



. ‘ < 
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r. 



r, 



13. Vectors u = £l, 2 , 3 J and v = 



1 
2 



3 



-fci t 



.>> . 

T ■% 



3^ © • • • • "tV ••••••••••••••••• 

(equal/unequal) 



* $ 
£ * * * 

% 9 * 

- ■Vjf v. 



I 



£ 



I 

l 

I 

| 

S' 

I 












1 



14. 



, * # A *' 
^ ; v 



ifc. 



& 






£ 



unequal 



?a » ? 



Let u « Fu^, Up, U3I and v = I v^, vp, V 3 I be^two row vectors, 

*" “ ^ #' 5 

% 

Then their sum u + v will be component -wise addition as follows: 

l l + v i> u 2 * v 2> u 3 + 
and v o £ 4 , 5, 6^ 



u + v « [u] 

B»» f 1 ’ 2 ' 3 ] 



& , j • 



a.. , * 
w> ; 



Then u + v =* J , 






1 :.*;■* 

a 1 % * 








' J 



M „ 
\ 



44 Vi 
**■ 



9 * 

» > 



j/ ' 

#v* 

•a 4 - ** 
' K \ 

* * 



[5, 7, 9] 



15. Let 





r 

u i 




" v i " 


U 0 


u 2 


and v *= 


v 2 




> 

. u 3 J 


i 1 


. v 3„ 



^r. 

f . ** % 

4 ' * }& ^ 

jfe * £ / 
*' 

4 ji. i ' 

t V. ’ .* % 



Then u + v 



u i + v i 
<»2 + V 2 
U3 + V3 



* k - ; ■ 

* ’-*<■ ' 

/. 

* i* /’ » 

jjiSW; " 

is a vector of siz# -; L4 






ig** 1 



T£* m< 

* i\ ,> 

' *&\V 
*» •■ 



I 







Then u + v » 



• • • • • 



1 

5 

9 



17. 



Let u e £ -2, 5, 2jj and v = £ 3 , -1, 4 J 

Then u + v = ^ , ~| 

and v + u a £....., ...... J 

and therefore u + v and v + u are •••••••••• 

(equal/uneaual) 



[l, 4, 6 

[l, 4, 

equal 





* • • . •> .,4^^, „ '_ „ „. »«....»■* 



*1 
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18. Let 


f* nr> 

3 




» 

-2 


u ** 


-2 


and v as 


3 








-4 

"» w 



Then u + v » 



and v + u ■ 



and therefore u + v and v * u are 

(equal/unequal) 




19* The sum of two row vectors of the same size 3 is the row vector of size 



3 






20* The sum of two column vectors of the same size 4 is the 
of size •••••••••• 



vector 









I 



I 



l 



I 



% 



I 



i 



$ 

% 



t 



i 



I 



s <• 



22 . 
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column. 



21, The sum of two n-component vectors (either row or column) is obtained by com- 
ponent-wise addition, and yields another n-component 



vector 



COMMENT 



Observe that we do not add vectors unless they are both row or both col- 
umn vectors, having the same number of components. 

Also we have seen that the order in which vectors are added does not 
matter; that is, 

u + v = v + u 

where u and v are both row or both column vectors. 

Therefore commutative law of addition holds true in the case of 
vectors also. 



The multiplication of a number a times a vector v is defined by component -wise 
multiplication of a times the components of v, For the three-component row 
vector we have 



av - a [v x , v 2 , v 3 ] • [ av x , av 2> av 3j 

and for the three-component column vector we have 



au « a 



U 1 




• • • • 0 


u 2 














L J 



I 



> 







i 



r 



i 



i 



K 



8 



I 



i 



% 



i 






23. If v s Jl, -4, 7, 9 J then 

0 [ » M, «; •••••, | 


! 


[ 5 > 


- 20 , 35 , 45 | 


24. If u = £ u i> u 2 , U3 J , then 




(“l)u = (-1) [u x , u 2 , u 3 J 




a r. * 


1 


| ^ | 


[-U1, 


, -u 2 , -U 3 ^ 


25. If u is any vector then its negative -u is the vector (-l)u. Therefore if 


u a lux, u 2 , u 3 , u 4 J , then 




. r -» r 




-u= (-1) Lu x , U 2 , u 3 , u 4 j = l 




[- 


u l> * u 2> " u 3» -U4J 



£ 



then 



26. Let « 53 |"s, 6, 7J and v = £ 1, 2, 3*] , 

«- [». «■ • f». «, 3i ■ [. ...... ] 



eric: 



m 



Mi 






m§ 









aaMMaafe 



”.»r . ~ " ; ’. "• ' "' • -:" w ; rs /":•" r -- •, ^.v: --•, .•>•;' ■ :. •';. .'. ..’" : ' 
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27. 



Let 





v 




*»■ j. 

V 1 


U «S 


U 2 


and v = 


V£ 




- U 3- 




v 3 

» !»» 



then 



i 


'»xl 




» 

< 

l-» 

i 




- — 
•••••••••••• 


U-V SS 


CM 

9 


- 


CM 

> 


85 






U 3_ 




v 3 

** m 







Uj^-Vj 

u 2“ v 2 

u 3“ v 3 



Let 


* 

*-* 

O 




“-3 " 




■' a. 

1 


U O 


8 


> v * 


5 


, and w « 


3 




- 6 „ 


! 


7 ^ 




.4, 





~10~ 




-3 * 




M «■ 

1 


u+3v-5w ts 


8 


+3 


5 


-5 


3 




. 6.. 




7 




*4- 



Then 



Tuf 




p* ** 




gr •* 




f' »« 


8 


+ 




+ 


; 


85 




.6. 








> 

*■**' * M 







o 

sic 



KBiiiltiiHHBiBSlitaMilBMHfii 






**^**** “ ~ — ;■■ ~ ~ — - - . . ~X, \U-~.i- J .- ,^. ■ — - -^ 
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» 



r -» 


, 


• -*■ 




-9 




-5 


1 l 

| 


15 


> 


-15 


l 


21 

-*• « 




^-20 _ 


I 



4 

8 



7 



r 




-i r 




1 


r «n 




29, Let u « 1 


7, 0, 


-5 J , v - [ 2, 


-3, 1 


j , and w =* 0 , -1J , 


Then 


2u-3v+5w 


“ 2 [ 


7, 0, -5] -3 


[ 2 , -3, 


1 j + 


5 fi, 0 , -ii 




f* 




T r 


U 




r 


*m 


83 I • • • • } 


!••• ) 


••••I ^ !••••; 


• ••• ; 


• • • • 1 


H* I •• •• ; •••• ; 


• • • • 


=[...., 


• Ml; 


....] 














[ 14 , 


0 , - 10 ] 


+ [-6, 


9 , - 3 j + [ 5 , 0 , 


-] 



[l3, 9, -IS] 



30, Is the sum 




possible 



(yes/no) 











:/ 

% 

e 

I 
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31. A vector whose components are zero is the zero vector. 



The vectors 



0 = 



0 

0 

0 



and 



0- [o, 0 , o] 



are three-component 




5' 



i 



| 



I 

I 




$ 



i 

I 

I 

I 

k 

f- 

l 



zero vectors 



32. Let 




, then 




33. 




Let v 53 fv^. 


v 2» v 3] 


v + 0 « 


| v l» v 2> 



9 then 

v s] + [°» 



0, 



0 



] 



[ 



# • # • ) 



• • • • 5 



• • • • 







I 



fe. 



I. 



% 



% 



t: 






t 



i 






& 



? 



w 
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38 



. Let u « J 2 , 6, l] and v « 



3 

4 
7 



f-3 



Then u v 



r 



~-l 



2 9 6, 1 



>1 



17 



• # • « 



6+24+7 =37 



39. 



Let u * [ u l> u 2* u 3 j an< * — 



V 1 

Vo 



/3j 



Then uv ■ jjij, u 2 , u^j 



V 1 

v 2 



L 3J 



u^v i+U2 v 2^ u 3 v 3 



40. We figured out in the previous frames that the product of a row vector and 
a column vector (both of the same size) is a 



( number/ vec t or/mat r ix) 



number 



<5 

ERIC 



MB&SBSlBfii 



iAi 



yiffliSiicMBaaaMii 






* '* 4 ' " *^ '*** ’^■v*?.-" *~ ' ~ ’ J ^“ ' *~— -~ 

- W 1 t & "'-' , ., 






**s*+-\ 






*$t ,# j* « *’ 

'#* 3f ■ # r ■ 

V*' “**••• 

,-v: 

, ■$■ 4; 4 , **■- > ■ si*V JVV 4 'A 2, 3& 

8&&9S'' %* * 

’Sn ^'7; * ^ 

v ‘ SiSr 
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*r 






#7- ‘^t^y**-*-* 
k, i^;»*i*" S! •■*' 
t* 

3* piyftM . . 



4* W 



* ••«&- * v*" 
>#1" ... *■ ; ,' 

■ ^v j V£ ; * 

,>!i : ** -.«y tt 



‘- * 
. ^’ r 



% 4ft 



fc?» 



v ^ 

siE*. '-•* 



**£"■'*&'** *** 



,'V' > 



fcdBr* ■ »j$v**- * 



*$§** , , W ,. 



. <Pfci “ 






At * 

u 



» W - -4;.#, -*■ 

-V'S-jt - -t-A , *, .v * 

*ft&: 4 # **#* . * 

Jk* •„*:f. * •;*«*' * #„ 

‘3*%i\**'* $ < '+ ' ' ' .- v , 

~ > ’ ■ "*? ** -■ ■'•* % 

> f * ** .1 

*Y ?£|*‘ «* *. / 

4 ., A . ■ r € .*> ,* w ft 

:■.. sp* » ^ ” <4 a* 

f* ‘.\V # ‘- >, - , 

** 

"I. * • y 

% * ? 



fv* 1 -m 

i , r -h \ 

n 

, ^ 




v l” 


u 0 

_ 9 »2 


and v ■ 


v 2 


; u 3j 




v 3 

-• J 



V. *f 
n A* rh ; , 



Then 'u v 

‘«M 



# 



ii * 



— • 
U 1 




v r 


u 2 




v 2 


. U 3- 




_ v 3 — 



- 4 -» 

- ■ «** ^ 

~ i * 

^ 

k * ,7 , : . : ■ 

_ 

•« V # ?< 
•% ^ 

•vt 4 #. 

* ■ <" -f ,) 



& 



* > , 



> ^ ^ ^ ** 
’; ^.sj^ ."’' ■’ 

> ./> ‘ 4(. 

V .-• » * . 

■r&7> 



are not 



11 ■*'*, 4- , 

m*h', $ 

. k* ' ■ ■*. 



,> 



, 1% 

* . a 

^ * # * 



*• 



Coes it make any sense? 



(yes/no) 



'V* - 4 '*; 

* . 7 * 



^^5 * 



-. ■* *• ^ 

, >• * St* 

\5 % 

£<r => 

«. ' •• 

xv - 

. :«%> 

irtU 



no 



* .« ■ ^^• , 

4 ^ *4'* 4 i< 
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46* The multiplication of a row vector by a column vector is called a minor 



product or a scalar product* 
The following product 



[v u 2> u a] 



v i 

V 2 



L V 3J 



1 eiffJ £* 



.A 

A ^ / 

= # 5f ' 

( ■ 4 

W’'? 



(is/is not) 



a minor product* 






^ #'* * 4t 

• ’ Jf > 

*^v- - ; 

7fVt * 

4 ^ / * 

" % -:• a3f 2 

i *s 

■ w T&*ji*'" 

■ . #•% ; 



is 



’> tfe 



^ O : 

ERIC 






ISI 



I iiij 

n: 

H 

I 

Jl 



(f 



ft 



f 

i 



a 



. . . .«w mct s ^ t tt&^.VUtt^VP-X.'' 
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47. We learned that a matrix can be transposed by changing rows into columns and 
columns into rows. The same principle can be used for transposing row and 
column vectors. Therefore after transposing a row vector it will become a 



vector and a column vector will become a vector 



column 

row 




48. 



Let u and v be twc column vectors of equal size, then the multiplication of 

the transposed of vector u by the vector v conformable for 

** (is/is not) 



multiplication. 



is 



49. Let u « [«,, U 2 > u^, u^J and v «* jv^, V£, v^/ v^J . Also v* is the 

transposed vector of v. Then 











50. 



51. 



ERIC 






<216 



Let u » ju^, 
transposed of u. Then 



V v 



u,l and v 

4 j 



V- [n» V v 3» v 4] 



AI30 u is the 



r 



S!" [ v l> v 2> v 3> v 4*j 



“1 
U« 



3 

u 4 



• • • 



• • 4 * ••••• * 1 * ••••• * 1 * 



• •••• 



v 1 u 1 +v 2 u 2 +v 3 u 3 +v 4 ,: 4 



I 



$ 



I 



^1 



'5 



•: 



/ 



i 



% 



i 






Let u 



— " [v V "a] 



, E ■ 



V 1 

v 2 



/ / 

and u, v their transposes. 



r 3 1 

J j 



Then 



EV - [v «2> u a] 



V 1 

v 2 



rv 



andvo.^i, v 2 , v 3 J 



'1 

“2 



Ai 



' / 



Therefore u v and v u ••••••••••••••••••• equal. 

(are/are not) 



























! 



IB 








m 



‘ ’ *- ■* •...-’• ' .— ..--pr^ r " ' 

~ -*—• — » — ■■ ■■ • ■ ■ — . — „ — — „ 



'i- 



tt 



£ 



I 



u 



& 
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u 1 v 1 +u 2 v 2 +u 3 y 3 



v l u l +v 2 u 2 +v 3 U 3 
are 



It is true In general that the minor product of two vectors is equal to the 
product of their transposes in reverse order* 



52, The second kind of product of two vectors is called the flavor product "of two 



vectors. The major product of two vectors is the product of a column vector 
by a row vector^ . The foiloWipg product 

i**” 

' u. 



u 



2 [> v 2» v 3» v 4] 



u 



3 | 



is a ,,,«•••••••»••••••••••••• product of two vectors, 

(minor/raajor) 



major 

Notice in major product two vectors need not to be of the same size. 



53, In a major product, two vectors of different sizes, the first a^ol^unn- vtctor 



and the second a row .vector,,, conformable for multiplica- 

(are/are not) 

tion. 



are 




1 



$ 

i 



I 



% 

% 

i 



4 

1 



% 



i 



3 














$ 



| 

I 



I 



l 






- *J, .^. 5 ^M.4U!IWr ’ ' J 



tauacsasa^^ 
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54. The major product of two vectors can be performed in the same way as the pro* 
duct of two matrices, i. e # , multiplying row by column. The major product of 
the following two vectors will be 



3 

1 

5 



[l, 4, 2, ft] 



© • • • # 



12 



18 



6 



20 10 30 



55. Multiply the following two vectors 



Ui 



u 



2 

u 3 



[v V v 3 ] " 



“aJ 



I 



l 



% 




u l v l 


U 1 V 2 


V 3 


U 2 V 1 


u 2 v 2 


V 3 


U 3 V 1 


V2 


U 3 V 3 


U 4 V 1 


V 2 


V, 







J '“ ' - ‘ I " ■ ■ » ', * “‘ ■ • ' *S?-T 1 ■-. '.. . ■■»,-. ■»-'• . -'■ — : . .. f | J | r . ’ , ' * '* 'll' *' ' " : ' • ' ' * •" r .'!-' ' ' " . _:’l. ’ "' ■?■' li 






56 
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• Let v ■ £l, 2, 3^ • The minor product of v and its transpose v is 



s _ r 1 

V V s j • • • • • f 

(write two vectors) 



2 2 2 
1+2+3 = 1* 



3j 



2 

3 



57# Let v « V£, v^"] . The minor product of v and its transpose v is 



V V = JVj, V 2 , VjJ 



V 



V, 



L v 3. 



V 1 +v 2 +v 3 



56. The minor product of a vector by its transpose is simply the sum of 
of the elements of the vector. 



squares 



59, Let v * b> 2 > 3 1 • The major product of v and its transpose v will be 



/ 

v v 



1* « 




2 ”1 

l z 2 3 




r i = 


VO 

CM 

CM 

CM 




L J 


, 2 


rf MM 




3 6 3 J 



(write two vectors) 



s 









i 









I 






n 



I 






I 



I 






I o 
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1 

2 

3 



2, 3J 



60. 



Let v o £v^, v 2^ v jJ • The major product of v and its transpose v will be 



/ 

v v 



V 1 

v 2 


[ v l> v 2> v 3] 


a 


MM 






■ v 3. 








• 9 • t • • • 





V 1 V 2 



v l v 3 



V 2 V 1 



V 3 V 1 



V 3 V 2 



V 2 V 3 

2 






COMMENT 



A square matrix in which the corresponding elements of corresponding rows 
and columns are equal is called a symmetric matrix. As we saw in the previous 



frame that the major product of a row vector v and its transpose v^ 



/ 

v v 



— * 




2 I 


V 1 




V 1 v l v 2 ^3 




r i 


2 


v 2 


\?l> v 2 , V 3 J * 


V 2 V 1 v 2 V 2 V 3 


73J 




_ V 3 V 1 V 3 V 2 v 3 2 _ 



has the elements of the first row equal to the corresponding elements of the 
first column, elements of the 2nd row equal to the corresponding elements of 
the 2nd column and the elements of the 3rd row equal to the corresponding 
elements of the 3rd column. /lit- general the major product of a vector and 
its transpose is a symmetric matrix. 



’ 1 



8 



6, 

■i 



a 

'I 



I 

£ 

5 






$ 

I 



% 















I 



i. 
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61. Let 



u 



u 



1 
u 2 
Ua 



. Then the major product of u and its transpose is 



/ 

u u = 



»r 



u, 



u„ 



[l- u 2» u s] 



and the resultant matrix is a .. * matrix. 



(kind of matrix) 



2 



U 1 


u l u 2 


u l u 3 


U 2 U 1 


2 

u 2 


u 2 u 3 


U 3 U 1 


u 3 u 2 


«3 2 




symmetric 





62. The vector being a special case of the matrix, the commutative law of multi- 
plication which does not hold for the product of two matrices, also does not 
hold for the product of two vectors* The minor product of two vectors is 



a *•••••• while the major product of the two 

(number/vec tor/mat r ix) 

vectors is a It is true because the 

(number/ vector/ matr ix) 

commutative law of multiplication does not hold for the product of two 



vectors. 



number 

matrix 



3 



% 



J 



■3S 



5 



I 
















,W^w~,W,U, 



w-aiwiaowas*^^ Beswajoc; ssgaaaswtsawi 
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2. Lee 



x + 2y + 3z ■ 1 
4x + 5y + 6z *= 3 
7x + 8y + z ** 5 

The above linear equations can be represented in the form of Au = v as follows 



••••• ••••• ••••• 




’X* 


i 


fr 






y 


= 




m * • • • • #•••# ••••• 




z 

td 




L 9 



1 

4 

7 



2 3 



5 

8 



6 

1 



y 

z 



1 

3 

5 



3® Represent the following system of linear equations 



a^x + b^y +> CjZ 



a 2 X + b 2 y + C 2 Z 



= d 



a 3 X + b 3 y + C 3 2 



* cl. 



in the form Au = v 



“ ^ 








** «a 






• ••ft 


a 


► 


w 








*- a 



"*1 b l *{ 




x* 




V 


a 2 b 2 c 2 




y 


SS 


d 2 


a 3 b 3 c 3 




z 

« 




, d 3. 



ERJC 



..MifllifflMlTlllilJ 






| 

I 



|? 

| 



i 

I 



a 

] 












i 











4 . 



5 . 



6 . 
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We have seen that a system of linear equations can be represented in the 
form Au « v, where A is the coefficient matrix for the system of equations, 
and u and y are vectors. Also A and u are conformable for multiplication 
and the product Au is a 



(vector/matrix) 



vector 



-1 



If A is the inverse of the matrix A, then the equation Au ■ v can be 
solved for the unknown vector u, getting 



jl . 



A -1 v 



Therefore a system of n .linear equations for n unknowns may be written in 
terms of its coefficient matrix A as the equation Au « v for the un k nown 
vector u. If A has an inverse matrix, then the solution is given by 



u 



-hr 






:«v. 



i 
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9. 



If 



x 

y 

z 



-2 

2 



Then the values of x, y and z are 



x 



y ■ 



ano z = 



3 

-2 
2 



3 , " 2 , 



5 ! ' 

! ; 

I 

i . 



10, Let 



4x + 5z » 



y - 6z = -2 



3x + 4z * 



The above system of linear equations can be represented in the form Aufv as 
follows 
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Unit XVI: The Theory of Linear Dependence 



1 . 



We learned in a previous unit that a determinant is zero if any of its 
rows (or columna) are linear combinations of other rows (or col umn s) , The 
determinant of the matrix 




is zero because the column 2 is equal to times of column 1, 



3 

(three) 



2 • M in the matrix of the previous frame we denote the columns 1, 2 and 3 as 

vectors v^, v^, and v^ respectively^ then: 

(a) .^2 “ k -1 w ^ en k 0 ••••• 




& 



j 

i 

s 



.3 






I 

1 







3 



(b) but Vg £ 



v, whatever the value of k. 







then they are called linearly dependent and if they cannot then they are 



In the previous frame linear dependence and independence has been defined 
in terms of only two vectors. The general definition applies to any number 

of vectors. Thus in general, n vectors all of the same order, v^, v 2 * , 

v n are said to be linearly dependent if non-zero constants kj^ k 2 , 

k n exist such that 



If zeros are the only values of the k's for which this equation is true, the 
vectors are said to be linearly independent. 

3. In the matrix 
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(c) When two vectors of the same size can be represented by a linear equation 



called linearly independent. Therefore, in the previous frame v^ and v^ are 
and v^ and are 



linearly dependent 



linearly independent 



COMMENT 



k lD- + k 2— 2 + 




13 2 



A a 2 6 7 



5 15 11 



the columns 1 and 2 are linearly 
A’ 1 does not exist. 



, and | A| a 



, therefore 



am 



- * 



wUi : v?* ! „ >> * ; ' '» w <*».'; V- » *. > Ir'^, " '^* t>i ^7 > • w £*■'• • ” > si ? W , v?;i »*< i J. : -< ^ J^.Vi J-A 1 ut> ^ JvsWr-l .^-:jW'-=' AW r v • 
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dependent 

zero 

4. In the matrix 

13 4 

2 7-1 

5 9 11 

the columns 1, 2 and 3 are and also |A| £ ...., 

therefore A”* does exist. 




1 



I 



I 



*1 



I 

l 

I 



I 



f 




■l 

tr 






3 



•r 



| 
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linearly independent 
zero 



COMMENT 

The property of the determinant which we noticed in the last two frames 
for a 3 x 3 matrix, is true in general: that if the rows (or columns) of 

a matrix are linearly dependent its determinant is zero and its inverse does 
not exist; conversely, if they are independent its determinant is non-zero 
and the inverse does exist. The minimum condition for the linear dependence 

of rows (or columns) of a matrix is that one of them be a linear combination 
of the others. 
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5. The number of independent rows in a matrix equals the number of independent 
columns and vice versa. The matrix 



-1 



5 9 11 

has 3 linearly independent rows therefore it has 
pendent columns also. 



linearly inde- 



6. The rank of any matrix is the number of linearly independent rows (or 
columns) therein. The matrix 



A 



7 -1 



9 11 



has 3 independent rows (as well as columns). Therefore the rank of the 



matrix A is 










1 



* 



& 



I 
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COMMENT 



Let A be m x n matrix. If A contains at least one r-rowed minor deter- 



minant that does not vanish, but no non-vanishing (r+1) -rowed minor deter- 
minant, A is said to be of rank r. If A * 0, the rank is said to be zero. 
If r * m or r » n, A contains no (r + 1) -rowed minor determinant. If r<m 
and r<n, the definition implies that A contains at least one (r+1) -rowed 
minor determinant but that each such determinant vanishes. 



7. From the previous comment it is clear that the rank of the matrix A of size 
m x n where ra<n is equal to or less than 



m 



8, If A is a square matrix of order n and its rank r < n, then | A| » 0, and 
A \ exist, • Therefore, the matrix A is singular. 



(does or does not) 



does not 



9, If A is a square matrix of order n and its rank r=n, then |a| £ 0, and A~* 



(does or does not) 
and is called matrix of full rank. 



exist. Therefore the matrix A is 



(singular or non-singular) 















234 



does 

non-singular 

Note: The rank of the inverse of a matrix A of order n is also n, and 

I A’ 1 ! =1/|A|. 



If A is a matrix of order mxn and has rank r 9 then A has r linearly 

••••••• rows (columns). 

(dependent or independent) 



independent 




I 



1 



I 

i. 



i 

I 






% 

l 

I 

S 

?{> 



i 

I 

I 



11* To find ojit the rank of the matrix 



A 



13 4 

2 7-1 

5 9 11 



(a) the first row is multiplied by -2 and -5, and then added to the 2nd 
and 3rd rows respectively to make the sub -diagonal elements of rows 2 and 
3 in column 1, equal to zero; we get 
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I 
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(b) The second row of A^ is multiplied by 6 and then added to the 3rd row 

of A l to make the sub-diagonal element of row 3 in column 2, equal to zero; 
we get 



A/ 



-9 



4 



0 



-63 



(c) In &2 we have made all the sub-diagonal elements of A equal to zero. 

Now the number of non-zero diagonal elements of A 2 is the rank of the matrix 
A, Therefore the rank of A is . 



12. To find out the rank of the matrix 



A = 



8 



-1 



(a) the first row is multiplied by .... and and then added to the 2nd 

and 3rd row respectively to make the sub-diagonal elements of rows 2 and 3 in 
column 1, equal to zero; \<re get 

(continued on next page) 
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COMMENT 



The method of finding the rank of a matrix, which we have used in the 
previous frames is quite general: initially the first row is used to reduce 

the sub-diagonal elements of the first column to zero, basing the calculations 
on the first term of the first row. Then the elements below the diagonal 
in the second column are reduced to zeros using the second row, based on its 
diagonal term. Because the sub-diagonal elements of the first column have 
already been made zero in the first step of this procedure this second step 
does not effect their values. And so the process is continued, using each 
row this way in turn, until all remaining rows are zero or until the last 
row is reached. The number of non-zero diagonal elements is then the rank. 



IS 



SMB — MmtmMmsmmm 



j 





iiiin i >■<<— a in in rin ii Mm I i 



Ills" ^ ^ * 1 ^ Wi ^ L...!.>, .L-,. I 






' C* 



*" V. 



, ❖ 



, * 



233 



Unit XVII: The Characteristic Equation of a Matrix 



1. Let 



A 




* 



¥ 

,# 
f# 

*4? 

. V 



£ 



A-kl « 

* *> « 

" * Yi ‘ ^ 




** 



be a square matrix of order 2 and k is a si|lar variable then 

0 

I 



if 



/. * 



* * ' 



% 



* *4 



. '• 



’* *\’*£ 
4 > 



* V#* 

'+% 



>** 



i«h* 



' *#* 



A „ a i 

■ *%T ? 



. * 



** * 



■# **„ 



V,t: * 

> * 



» . % ;; 

- . •* +i‘& 



3-k 

I 



5 

7-k 



* r \ *; * 



v* " * 

# arin -square matrix and A a scalar variable, then 

* * - * - *#> ■;*** * *** » ■ * F 
, ¥** 



# 



* * v * ' 

'■ * ,. 



( a ) matrix B-A I is called the characteristic matrix of B. Therefore 



4 

^ if 



« * m * r _- 

,* C* v„S 5 



Ip' the .previous frame the characteristic matrix of A is 



* „ ■ $* 
% * J* ■■ ■ 



> *». l 
* 



* «* 

£ :t*< ;•'* * 
* ■,;* , 



** 

, * i 

*A & 



f 3 
as* ^ 






* f f 



t- 



‘ .^c A 

aUK V*; 



•* , , 

* 



3-k 5 

1 7-k 



'If 

•1*.v. 



V f 

Sv 



* * f ^ 



* 



. A 
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a 
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I 



I 



; 






j 



1 



a 






? 














(b) the determinant |B-Al| - f(A) is called the characteristic determinant 
of B; the equation f(A) ° 0 is called the characteristic equation of B a 
From frame 1 we have 



| A-klj = 



3-k 5 

1 7-k 



(3-k) (7-k) - 5 • k - 10k + 16, 



therefore k - 10k + 16 = 0 is the of A. 



characteristic equation 



(c) the roots of f (A) s 0 are called the characteristic roots . the latent 
roots » the characteristic values , or the eigenvalues of B. Solving the 
characteristic equation of A, k^-10k+16 =* 0, k » 8 or 2« Therefore the 
characteristic roots of A are «••••••• and 



8 and 2 



3 . 



Let 



A 



14 1 

2 10 
-13 1 



and i 



(a) The characteristic matriK of A is 




0 



0 

I 0 




f 



0 ! 

\ 

i 




0 
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4. Let 




4 -2 

2 -2 
-2 -1 



(a) The characteristic equation of A is 



| A- M| = 



0 



2-A 4 -2 

4 2-> -2 

-2 -2 - 1 - > 



0 



(b) The characteristic equation of A is 
2- A 4 -2 

4 2 -A -2 

-2 -2 1-A 

which after expanding the determinant, becomes 
-,X 3 + 3* 2 + 24> + 28 a 0 
or (>.+ 2) 2 <7 - ,*) - 0 



and 



2 



> 



2 



3 



7 



Therefore the characteristic roots of A are ...... 



• • • • • 
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5. Let 



-[::] 



and its latent roots are 2 and 8. Let u 



V 




V 




and v » 


X 


"2. 




- V 2- 



be two vectors 



such that 

Au = 2u 



and 



Av = 8v 



or 



E 313 • 'Cl - E 3[3- •[::) 



or 



t:.3 • M - 1:::::::::::::] • [:;] 



h* s -»l.M ^ p'l ♦»’«]. M 

L U 1 + 7u 2 J L 2u 2 J L v 1 + 7 v 2J L 8v 2j 



From the previous frame we have 
3 u l + 5uJ 



Uj^ + 7u 2 , 



ca 



and 



3v^ + 5 v 2 



V X + 7v 2 



Cl 



Taking the arbitrary values of u^ « -5 and * 1 ; and then solving the 
equations 3u^ + 5u^ ® 2u^ and 3v^ + 5v^ « 8vp we get 
u 2 a ••••••• and v 2 « ••••••• 



•• 1 






I 



I 



£ 



and 















f- - ' " -^ ^ ~- v^W. •sv^a,»^U^ *| 
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From the frames 5 and 6 we can see that corresponding to the latent roots 
2 and C of the matrix A ® | | , there are two vectors 



ri n* 

f-s] rn 

u a I and v » I 

UJ ' UJ 



such that 



i* ©• 



Au 

h 

3 

1 



and Av 



8v 



- 2 



and 



• 2u 

:] it 

Here u and v are the latent vectors of A corresponding to the latent roots 
(2 and 8) of A. 



to 



[: B • •[:] 



8# Let A be a matrix of order n and therefore the characteristic equation 

| A-Al| * 0 is a polynomial in A of degree n. If it has all the solutions 
then: 



(a) How many latent roots of A do you expect? 



n 



(b) How many latent vectors of A do you expect? 






4 









j 



| 






I 



1 



I 









4 









I 



E 



8 






i 



P 



EJ 
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(c) If ^ and u^ (for i 
of A, then 

Au^ = 



l, 2, ..., n) are the latent roots and vectors 






9, Let 



A 



-1 



-2 



0 



(a) The sum of the diagonal elements (which is called the trace) 
of A = ••••• 4- 4* = 5 



2 4- 3 4-0 



(b) The sum of the minors of order 2 of the diagonal elements (which are 
called the principal minors of order 2) 2 f 3 and 0 of |Aj 



■3 
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l 













1 



I 



£ 

me 









245 



2 -1 

3 3 



: 2 



(c) The characteristic equation of A is 

2 -> -1 1 



| A- M| * 



= 0 



3 3-> -2 

4 1 -v\ 

and it can be expanded by diagonal elements (called the diagonal expansion) , 



It is 



a v 2 

(- A) J +(-/v (trace of A)+(-A) (sum of the second order minors of A)+|A|» 0 



or 



-> 3 + A 2 (....) - >* ( 7 ) + 3 



(d) The roots of the characteristic equation 
-A 3 + 5^ 2 -7^+3 = 0 
or (A-l) 2 (>v-3) = 0 



are 









, and ^ 



l, 1, and 3 



COMMENT 



The method of diagonal expansion of the characteristic equation in the 
previous frame is quite general and can be used for a matrix of order n # 
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Unit XVIII: Special Types of Matrices 



1. Let 







a ll 


a 12 


a l3 


II 

< 


a 12 


a 22 


a 23 






a 13 


a 23 


a 33 


(a) 


The elements of 


the 


and 


the 


elements of 


the : 



9 • 1 



J • • • f 



» • • • • J 



*; • 



a ll> 


a 12> a 13 


a ll> 


a 12> a 13 



(b) The elements of the second row are 

and the elements of the second column are 9 . ... 



a 

a 



(c) The elements of the 3rd row are 
and of the 3rd column ...,., , 



12* a 22^ 


a 23 


12* a 22> 


a 23 


; t ••• • 



a 13’ V a 33 
a 13> a 23> a 33 



r- . £*-*x/r, * »• •* 



I 



% 

'I 



9 



I 
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(d) The transpose of A is 



A 



(name of the matrix) 



11 


a 12 


a 13 


12 


a 22 


a 23 


13 


23 


a 33 



A 



2. Now let 



B 



0 

-a 

-a 





a 12 


a i3 


12 


0 


a 23 


13 


-a rt#> 

23 


0 



then 



(a) The elements of the 1st row are •••••, •••••> ••••• 
and the elements of the 1st column are •••••; •«*••• 



°, a 12 > 

°, “ a i2^ 



13 



-a 



13 



(b) The elements of the second row are •••••; •••••* ••••• 
and of the second column •••••; • •••• 






m 






BS^MiagiBSSM 
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12' 


o. 


a 23 


12' 


o. 


" a 23 



(c) The elements of the 3rd row are 

and of the 3rd column . . . . . ^ . . . « . ^ ••••• 



" a 13' ’ a 23* ° 

a 13' a 23' 0 



(d) The transpose of B is 




0 


CM 

H 

« 

1 


-*» 




~ 0 


a i2 


pi 

»-» 

U> 

1 




a 12 


0 


-a n . 

23 


as - 


* a l2 


0 


a 23 


• -B 


a 13 

Wm 


a„ 

23 


0 




~ a 13 

W 


" a 23 


0 

m 





COMMENT 

We have seen in frame 1 that the rows of A are equal to the corresponding 

✓ 

columns of A, and also A =* A , therefore A is a symmetric matrix. 

In frame 2, the rows of B are equal to minus times the corresponding 

s 

columns of B, and also B » -B, therefore B is a skew-symmetric matrix. 
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In general a matrix A is said to be symmetric if it is equal to its 



transpose, i.e., if A » A or a^ a a^.; (i,j * 1, 2, ,,,*., n), A is 



said to be skew-symmetric if it is equal to the negative of its transpose, 
✓ 

j ^ 8 t ®£j ~ “ a ji* a ii s (^; j 53 lj 2, ••«•••, n) # 



3# If A is such a matrix that its inverse equals its transpose, i, e, A 



-1 



✓ 

A, 



then A is said to be an orthogonal matrix. 



Also AA 



-1 



AA 



or 



AA 



-l. 



(multiplication of A and A ) 



Let A be an orthogonal matrix and A a latent root, then 

|A ->I| =0. * 



Multiplying throughout by A we get 



| AA -A IA| o 0 



or 



j «••••• 



-AA| 



Multiplying throughout by 1^ we get 

a 0 



B- 1 * I 



because A a a 



A 



% 



? 



g 

§ 



3 



>S 

5 



5 



■ri 






I 



or l A -Fr° 

It shows that if is a latent root of an orthogonal matrix then so is l/^\ # 



I 







5, Let A be an orthogonal matrix so that AA<:» J.„ Then 



| AA| « 1 1| ** +1 or -1 



-T 1 ** 



! o* 



Therefore the determinant of A is 



Hi 



OsT, * 






. y . __W _ 



v ;r -T- 



'? s 



// * 
V 



+1 or -1 



6* Let 



-6 -4 

4 2 



*■ i ' 
* 1 * 
V 4 



V- 4 

' # 



3 

-6 

-4 2 -1 

be a symmetric matrix x*ith its latent roots A « -1, -4, 11. 



■ P ! 

, ** 



Corresponding to the latent roots 9 the latent vectors are 



•** 



1 

2 



-2 

d 



2 

1 

2 



and 



2 

-2 

-1 



* « 



* /? . . 4 



* a » it 



i.e. p 



p * # 



which can be represented as columns of a matrix P 

12 2 * 

2 1 -2 

-2 2 -1 

(a) P can be normalized dividing elements^ *of ' each column by the square root 
of the sum of the squares of the elements in the column. Therefore normalizing 
column 1 of P, elements 1. 2 and -2 should be divided 



♦ - V 

f , ■- 

t ** * 

' '* 



v*/* 



by \li 2 + “ 2 2 



2+2 

normalizing column 2, elements of 2, 1 and 2 should be divided by 

\j» • • 1 



3; 



3; and normalizing column 3, elements 2, -2, -1 



should be divided by \f .... + .... + .... 






I 

ft 



I 



1 



i 






§ 



j 












.... .^, j ; J,,i>. > i, ;V JIJI.J, 
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n! 2 2 + l 2 + 2 2 



,2 2 2 

x(2 + 2+1* * 3 



(b) Normalized form of F is obtained by dividing columns 1, 2 and 3 by 3. 
Therefore normalized P is 









IP* — 

12 2 


* 








P « 




= 1/3 


2 1-2 








-2 2 -1 


€ 


H* »4Ml 


' 



* 

P P 



1/9 



1 2 “2 

2 12 
2 -2 -1 



1 2 

2 1 

-2 2 



2 

-2 

-1 



1/3 


2/3 


2/3 


2/3 


1/2 


-2/3 


-2/3 


2/3 


-1/3 


of P* 

■m 


to P* 


is 



J 



= 1/9 



• • • • 



• • • • 



• • • • 



• • • • 



• • • • 



• • • • 



Therefore P* is an orthogonal matrix. 



1/9 



9 

0 

0 



0 

9 

0 



0 

0 
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